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SUMMARY 


A.  VERTICAL  ARRAYS  FOR  TELESEISMIC  oIGNAL  EXTRACTION 


Deep-well  seismic  data  recorded  at  UBO  in  September- 
October  1964  arc  usc'd  to  estimate  the  capabilities  of  three  3 -element 
vertical  arrays  and  one  hypothetical  6-element  vertical  array.  It  is 
shown  that  two  deeply  buried  receivers  can  be  employed  to  estimate 
accurately  the  signal  observed  by  a  surface  instrument.  It  is  estimated 
that  the  maximum  signal -to-noise  ratio  improvement  obtainable  from  a 
single  vertical  array  at  UBO  is  approximately  6  to  8  db  but  that  a  vertical 
array  should  be  able  to  provide  greater  improvement  at  a  high-noise  site. 

It  appears  that  the  percentage  of  ambient  noise  due  to  P-wavcs  is  a  dom¬ 
inant  factor  contributing  to  the  performance  limit. 

The  vertical  array  data  from  UBO  consisted  of  simul¬ 
taneous  recordings  from  three  seismometers,  one  held  at  a  depth  of  200 
ft  and  the  other  two  placed  at  all  possible  depth  pairs  of  five  deeper 
depth  positions.  Correlations  from  these  recordings  were  combined  to 
estimate  the  correlation  matrix  for  a  hypothetical  6 -seismometer  vertical 
array.  This  built-up  correlation  matrix  was  not  positive  definite,  i.e., 
was  not,  in  fact,  a  possible  correlation  matrix.  All  attempts  at  generat¬ 
ing  a  6-channel  multichannel  filter  system  failed,  both  in  the  time  domain 
and  frequency  domain.  The  correlation  functions  are  in  good  visual  agree¬ 
ment  with  theoretical  correlation  functions  and  indicate  considerable  man¬ 
tle  P-wave  noise  at  UBO.  However,  to  insure  the  necessary  mathema¬ 
tical  interrelations  between  the  auto-  and  crosscorrclations  for  MCF 
design  for  an  array  of  seismometers  a  simultaneous  recording  from  the 
full  array  seems  to  be  a  practical  necessity. 

B.  DECONVOLUTION  AND  AUTOCORRELATION  STUDY 


The  application  of  deconvolution  and  autocorrelation  tech¬ 
niques  to  three  ensembles  of  earthquake  records  from  the  multichannel 
processor  at  CPO  is  being  studied.  Statistics  cf  discrepancies  between 
observed  and  USC&GS  reoorted  origin  times,  magnitudes  and  depths  of 
focus  are  examined. 

Ensemble  anti  subensemble  average  autocorrelations  are 
presented  tor  Ensemble  I  (73  Kurile  Islands  events).  Very  little  differ¬ 
ence  can  be  observed  among  the  subensembles  or  between  unweighted  and 
normalized  averages.  It  is  hypothesized  that  certain  features  in  the 
autocorrelations  are  due  to  reflections  from  the  Conrad  and  Mohorovicic 
discontinuities  under  CPO. 


in 


Ensembles  II  and  III  (60  and  119  events  of  worldwide  dis¬ 
tribution  are  described  and  the  results  of  individual -event  deconvolution 
presented.  Although  signal  waveforms  are  contracted  somewhat,  there  is 
no  obvious  improvement  in  detectability  of  later  phases.  The  nature  of 
the  ambient  noise  recorded  at  CPO  is  discussed  briefly. 

C.  PROBABILISTIC  PROCESSING 

A  description  of  the  mathematical  operations  and  computer 
programs  required  to  perform  probabilistic  processing  is  given  in  Section 
III.  In  particular,  an  efficient  iterative  formula  for  inverting  a  multichan¬ 
nel  correlation  matrix  is  derived,  a  technique  for  generating  a  synthetic 
multichannel  time  series  with  a  given  correlation  matrix  is  discussed  and 
a  recursive  procedure  for  solving  the  generalized  eigenvalue  problem  for 
single -channel  correlation  matrices  is  presented.  The  status  of  these 
programs  and  future  goals  are  discussed. 

D.  STUDY  OF  PARTIAL  ARRAYS  AT  CPO 

Four  multichannel  filter  systems  for  CPO  were  designed 
using  the  full  array  and  three  partial  arrays.  These  systems  were  synthe 
sized  using  measured  noise  correlations  and  an  infinite  velocity  signal 
model.  The  MCF  systems  were  evaluated  on  the  computer  to  compare  the 
ability  of  the  partial  arrays  to  suppress  noise  while  preserving  high-velocitv 
signals. 


Results  show  that  the  partial  arrays  are  almost  as  effective 
as  the  full  array  near  1  cps.  Below  approximately  0.7  cps,  seismometer 
gain  inequalities  appear  to  be  cfiecting  the  output  of  the  systems.  At 
higher  frequencies,  the  K- space  aliasing  of  some  of  the  partial  arrays 
prevents  effective  noise  suppression. 

To  help  evaluate  the  effects  of  seismometer  gain  inequali¬ 
ties,  MCF  systems  will  be  designed  using  a  signal  model  incorporating 
statistical  gain  fluctuation.  This  will  be  the  subject  of  a  separate  report. 

E.  INITIAL  PROCESSING  AND  ANALYSIS  OF  THE  LONG  TFO 
NOISE  SAMPLE 


Section  V  contains  information  concerning  the  recording 
and  preliminary  processing  and  analysis  of  a  20-min  ambient  scLmic 
noise  sample  from  the  ring  array,  the  large  cross  array  and  a  3-component 
seismometer  at  Tonto  Fores*  Seismological  Observatory.  Also  presented 
is  an  analysis  of  calibration  data  for  the  large  cross  array  recorded  on  the 
same  day  as  the  long  noise  sample. 


F.  FURTHER  ANALYSIS  OF  TFO  AMBIENT  NOISE 


Section  VI  presents  the  complete  set  of  auto-  and  crosscorre¬ 
lation  functions  from  the  21  seismometers  in  the  large  cross  array  at  TFO 
as  found  from  the  long  noise  sample  discussed  in  Section  V.  Also  presented 
is  the  complete  set  of  10  coherences  between  the  two  horizontal  seismome¬ 
ters  and  three  vertical  seismometers.  The  analysis  of  this  data  shows  that 
1  he  4  to  6  sec  energy  is  coming  from  N60  E  with  a  velocity  of  about  3.4 
km/see.  The  other  energy  below  1.2  cps  is  mantle  P-wave  noise,  which  is 
definitely  nonisotropic.  Above  1.2  cps,  the  energy  is  random,  seismometer 
to  seismometer,  except  for  some  highly  coherent  lines  in  the  spectra.  The 
horizontal  seismometer  in  line  with  the  4  to  6  sec  energy  shows  a  high 
coherence  with  the  vertical  seismometers  in  that  frequency  range,  but  the 
transverse  seismometer  has  very  small  coherence  with  all  other 

seismometers.  This  is  further  proof  that  this  energy  is  Rayleigh  motion. 
The  presence  of  Love  wave  energy  on  the  transverse  seismometer  can  be 
inferred  from  its  lack  of  coherence  with  the  other  seismometers. 


vi 
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SECTION  I 


VERTICAL  ARRAYS  FOR  TELESEISMIC  SIGNAL  EXTRACTION 

A.  SUMMARY 


Deep-well  seismic  data  recorded  at  UBO  in  Scptcmber- 
October  1964  are  used  to  estimate  the  capabilities  of  three  3-element  ver¬ 
tical  arrays  and  one  hypothetical  6-element  vertical  array.  It  is  shown  that 
two  deeply  buried  receivers  can  be  employed  to  estimate  accurately  the 
signal  observed  by  a  surface  instrument.  It  is  estimate!  that  the  maximum 
signal-to-noise  ratio  improvement  obtainable  from  a  single  vertical  array  at 
UBO  is  approximately  6  to  8  db  but  that  a  vertical  array  should  be  able  to 
provide  greater  improvement  at  a  high-noise  site.  It  appears  that  the  per¬ 
centage  of  ambient  noise  due  to  P-waves  is  a  dominant  factor  contributing 
to  the  performance  limit. 

The  vertical  array  data  from  UBO  consisted  of  simultaneous 
recordings  from  three  seismometers,  one  held  at  a  depth  of  200  ft  and  the 
other  two  placed  at  all  possible  depth  pairs  of  five  deeper  depth  positions. 
Correlations  from  these  recordings  were  combined  to  estimate  the  correla¬ 
tion  matrix  for  a  hypothetical  6-seismometer  vertical  array.  This  built-up 
correlation  matrix  was  not  positive  definite,  i.c.,  was  not,  in  f  ct,  a 
possible  correlation  matrix.  All  attempts  at  generating  a  6-channel  multi¬ 
channel  filter  system  failed,  both  in  the  time  domain  and  frequency  domain. 
Tne  correlation  functions  are  in  good  visual  agreement  with  theoretical 
correlation  functions  and  indicate  considerable  mantle  P-wave  noise  at  UBO. 
However,  to  insure  the  necessary  mathematical  interrelations  between  the 
auto-  and  crosscorrelations  for  MCF  design  for  an  array  of  seismometers, 
a  simultaneous  recording  from  the  full  array  seems  to  be  a  practical  neces¬ 
sity. 


During  March  1965,  a  program  of  deep-well  recording  was 
conducted  at  Grapevine,  Texas.  This  report  presents  a  description  of  the 
data. 

B.  ANALYSIS  OF  UBO  VERTICAL  ARRAY  DAT/. 

1.  Description  of  Data 

As  described  in  a  previous  report  ,  a  program  of  digital 
recording  of  seismic  data  was  conducted  at  Uinta  Basin  Observatory  (URO) 
during  the  period  9  September  through  13  October  1964.  The  analysis 
reported  in  this  section  uses  data  obtaired  from  three  instruments  designated 
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SH-l,  DH-1  and  DH-2.  SH-1  remained  at  a  depth  of  200  ft  (0.  06  km) 
throughout  the  recording  period,  while  DH-1  and  DH-2  were  moved  about 
in  order  to  permit  simultaneous  recording  at  each  of  the  possible  pairs  of 
the  following  five  depths:  900  L  (0.27  km),  2900  ft  (0.88  km),  4900  ft 
(1.49  km),  6900  ft  (2.  10  km),  8900  ft  (2.71  km). 

Data  edited  from  the  field  records  consisted  of  50  records  of 
seismic  noise  (five  from  each  pair  of  depths  of  DH-1  and  DH-2)  and  16 
records  of  teleseismic  signal  (two  wrten  possible  from  each  depth  pair). 

Each  record  consisted  of  2  min  of  data  recorded  digitally  with  24-msec 
sampling  interval.  High-cut  filtering  was  applied  so  the  data  could  be 
resampled  at  72  msec,  and  gain  factors  were  applied  to  equalize  instru¬ 
ment  response  at  1 . 0  cps.  The  records  are  listed  in  Tables  1-1  and  1-2. 

2.  Preliminary  Investigation  of  Data 

Four  noise  samples  and  three  signals  from  two  sets  of 
recording  depths  were  selected  in  order  to  investigate  certain  properties 
of  signal  and  noise  as  observed  in  the  UBO  deep  well.  Playbacks  of  the 
original  records  are  reproduced  in  Figure  1-1.  No  visible  alteration  of 
the  records  was  produced  by  the  filtering  and  resampling  operation. 

Power  spectra  computed  for  the  selected  events  are  illus¬ 
trated  in  Figures  1-2  and  1-3.  Absolute  power  levels  n  iy  be  determined 
by  means  of  the  smooth  curve  which  corresponds  to  a  power  density  of 
2 

1  mp.  /cps.  The  average  system  noise  level  over  the  range  0  to  4.  0  cps 

,  ,  .  2  . 

was  found  to  be  -43  db  relative  to  1  m|J  /cps  at  1.0  cps. 

A  significant  difference  was  found  in  the  spectra  of  the  sig¬ 
nals  studied.  However,  only  relatively  small  differences  in  depth-dependence 
can  be  observed.  Some  depth-dependence  differences  would  be  expected  to 
result  from  differences  in  angle  of  emergence  of  the  signals.  The  peaks 
and  troughs  in  the  spectra  for  buried  seismometers  generally  are  located 
where  the  theory  predicts.  When  smoothing  is  applied,  the  noise  spectrum 
at  200  ft  is  found  to  be  relatively  constant.  Apparent  variations  in  the 
frequency  and  depth-dependence  of  the  noise  probably  are  not  significant, 
so  it  is  impossible  t  >  draw  conclusions  as  to  whether  important  changes 
may  have  occ’..  red  in  the  modal  structure  of  the  noise.  Nonstationarity 
of  the  signal  and  noise  is  the  most  probable  explanation  for  the  failure  of 
attempts  to  evaluate  a  6-element  vertical  array  using  the  UBO  data.  Refer 
to  subsection  B-4  for  further  details. 

For  each  of  the  selected  events,  a  short  deconvolution  fil¬ 
ter  was  designed  *o  whiten  the  spectrum  of  the  200-ft  trace.  This  filter 
was  apolied  to  each  of  the  three  traces,  and  coherence  functions  were 


1-2 


Table  1-1 


CATALOG  OF  UBO  DEEP -WELL  NOISE  SAMPLES 
(Dates  and  GCT  Start  Times) 


DH-1 

Depth 

(ft) 

DH-2 

Depth 

(ft) 

Noise 

Sample 

1 

Noise 

Sample 

2 

Noise 

Sample 

3 

Noise 

Sample 

4 

Noise 

Sample 

5 

8900 

900 

9/13 

9/14 

9/15 

9/15 

9/15 

21:17:20 

22:10: 10 

12:48:50 

14:01:50 

14:30:00 

89'D 

2900 

9/17 

12:07:50 

9/17 

13:10:40 

9/17 

1  3:52: 1 L 

9/17 

21:14:20 

9/17 

21:54:20 

8900 

4900 

9/20 

19:08:30 

9/21 

20:00:10 

9/21 
13:28:1  j 

9/21 

15:57:10 

9/21 

17:17:10 

6900 

900 

9/27 

18:55:20 

9/27 

20:47:10 

9/28 

14:56:20 

9/28 

15:41:20 

9/28 

16:42:30 

6900 

2900 

9/29 

22:16:20 

9/29 

22:54:20 

9/30 

22:07:00 

9/30 

22:45:26 

9/30 

23:01:10 

6900 

4900 

10/1 

19:17:16 

10/2 

20:28:50 

10/2 

22:47:20 

10/2 

23:34:10 

10/2 

23:52:50 

4900 

2900 

10/4 

15:34:10 

10/4 

16:05:10 

10/ 4 
18:43:10 

10/5 

02:08:10 

10/5 

02:36:30 

4900 

900 

10/6 

18:32:50 

10/6 

18:57:50 

10/6 

19:10:00 

10/6 

20:28:30 

10/6 

20:34: 10 

2900 

900 

10/8 

19:16:10 

10/9 
00:09: 10 

10/9 

00:23:50 

10/9 

14:27:10 

10/9 

15:00:10 

8900 

6900 

10/12 

18:58:10 

10/12 

19:  39: 10 

10/13 

02:14:10 

10/13 

03:53:10 

_ 

10/13 

04:09:10 

CATALOG  OF  UBO 


Ser. 

No. 

DH-1 

Depth 

(ft) 

DH-2 

Depth 

(ft) 

Date 

1 

8900 

900 

9/15 

3 

8900 

2900 

9/17 

5 

8900 

4900 

9/21 

6 

8900 

4900 

9/21 

8 

6900 

900 

9/28 

9 

6900 

2900 

9/30 

10 

6900 

2900 

9/30 

11 

6900 

4900 

10/2 

13 

4900 

290G 

10/5 

14 

4900 

2900 

10/5 

15 

4900 

900 

10/6 

16 

4900 

900 

10/6 

17 

2900 

900 

10/8 

18 

2900 

900 

10/9 

19 

8900 

6900 

10/13 

20 

3900 

6900 

10/13 

Table  1-2 


EP-WELL  TELESEISM  RECORDS 


Arrival 

Time 

Location 

M6 

A 

(degrees ) 

12:56:25 

Samoa 

5.3 

81 

13:35:54 

Unidentified 

Teles  e 

;ism 

13:37:07 

Peru-Brazil 

4,6 

58 

16:40:23 

Kurile  Is. 

4.8 

75 

16:38:05 

Honshu 

5.0 

82 

19:26:18 

Honshu 

4.8 

82 

20:33:51 

Kodiak  Is. 

5.  1 

32 

20:51:24 

Peru 

4.6 

55 

03:37:38 

Hokkaido 

4.7 

74 

0  3:46:45 

Hokkaido 

5.0 

74 

18:34:40 

Jan  Maven 

4.8 

57 

19:25:30 

N  Hebrides 

5.4 

94 

17:01:14 

Foxe  Is. 

4.  5 

40 

00:26:02 

Easter  Is. 

4.5 

75 

02:31:51 

Kurile  Is. 

5.  2 

71 

03:25:45 

Tonga  Is. 

4.7 

81 

UBO  VERTICAL  ARRAY  ORIGINAL  DATA 

24  MS  SAMPLE  RATE 
0  -  20.83  CPS 
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Figure  1-1.  Ambient  Noise  Samples  and  Teleseisms  Recorded 
in  the  UBO  Deep  Well 


POWER  DENSITY  in  db  RELATIVE  to  l(m±i2/cps)  at  Icps 


-20  h? 


FREQUENCY  (cps) 


FREQUENCY  (cps) 


Figure  1-2.  Signal  and  Noise  Power  Spectra  Computed  from  Data  Recorded 
at  Depths  of  200,  900  and  2900  Ft.  The  smooth  curve  corre¬ 
sponds  to  1  my^/cps 


POWER  DENSITY  In  db  RELATIVE  to  KmjiVcps)  at  Icpt 


computed  from  the  whitened  data.  These  results  are  presented  in  Figures 

1-4  and  1-5. 


The  coherence  of  the  noise  is  observed  to  be  dependent  upon 
receiver  separation.  When  the  separation  is  small  (200-ft  and  900-ft  depths), 
the  noise  is  highly  coherent  at  all  frequencies  to  ?..0  cps.  As  the  separation 
increases,  the  frequency  range  over  which  the  noise  is  significantly  pre¬ 
dictable  contracts.  In  the  case  of  the  largest  receiver  spacing,  there  is  no 
significant  coherence  above  0.  5  cps  except  for  a  few  peaks  which  do  not 
exhibit  time  -  stationarity.  The  persistent  trough  at  0.  5  cps  in  the  noise 
coherence  functions  for  6900  and  8900  ft  is  in  close  agreement  with  the  theo¬ 
retical  predictions ^  for  vertically  traveling  P-waves  and  with  the  trough 
observed  in  the  experimental  signal  coherences.  This  result  suggests  that 
mantle  P-waves  form  a  major  component  of  the  ambient  noise  at  that  fre¬ 
quency.  The  low  coherence  at  higher  frequencies  could  result  from  a  mul¬ 
tiplicity  of  modes  of  propagation. 

The  signals  are  observed  to  be  highly  coherent  over  the  fre¬ 
quency  range  in  which  there  is  appreciable  signal  power.  This  range  is 
seen  to  be  different  for  different  signals.  Notches  in  the  coherence  func¬ 
tions  occur  at  each  frequency  corresponding  to  a  signal  null  at  either  of  the 
two  receiver  depths.  These  troughs  are  found  to  be  in  good  agreement  with 
the  theoretical  predictions. 

3.  Investigations  of  3-Element  Arrays 

Experimental  estimates  of  the  performance  of  3 -element 
vertical  arrays  were  obtained  for  the  following  recording  configurations: 

•  200,  2900  and  4900  ft 

•  200,  6900  and  8900  ft 

•  200,  4900  and  8900  ft 

In  each  case,  the  five  noise  samples  and  two  teleseisms  were  used  to  deter¬ 
mine  the  (signal  +  noise)  statistics  as  shown  in  Figures  1-6  through  1-8. 

The  second  telescism  of  the  set  was  used  to  describe  the  signal  statistics 
for  the  right  side  of  the  filter  design  equation.  After  filtering  and  resampl¬ 
ing  at  144  msec,  optimum  multichannel  filters  to  enhance  signal  at  the 
200-ft  seismometer  were  designed  in  the  time  domain  using  both  the  Wiener 
least-mean- square -error  formulation  and  the  prediction-error  approach. 
Filter  lengths  were  25  points  (3.  5  sec). 

Evaluations  of  the  filters  were  obtained  by  applying  each 
multichannel  filter  to  noise  sample  1  and  telescism  2  of  the  corresponding 
data  set.  Power  spectra  were  computed  for  the  filter  output  and  trace  1 
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Figure  1-4.  Coherence  Functions  Computed  .iom  Data  Recorded 
at  Depths  of  200,  900  and  2900  Ft 
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5.  Coherence  Functions  Computed  From  Data  Recorded 
at  Depths  of  200,  6900  and  8900  Ft 


Figure  1-6.  3 -Channel  UBO  Data  from  200,  2900  and  4900  Ft 


-Channel  UBO  Data  from  200,  6900  and  8900  Ft 
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Figure  1-8.  3-Channel  UBO  Data  from  200,  4900  and  3900  Ft 


(200-ft  depth)  input.  Signal-to-noise  ratio  improvement  estimate'  computed 
in  this  way  are  shown  in  Figures  1-9  through  I-  11. 


These  results  ere  difficult  to  interpret  in  detail.  They  sug¬ 
gest  that,  in  the  frequency  range  below  3.0  cps,  a  3-elerr>ent  vertical  array 
should  be  capable  of  providing  an  average  improvement  in  signal-to-noise 
ratio  of  approximately  4  to  5  db  relative  to  a  single  seismometei  at  200-ft 
depth.  Since  the  noise  power  at  200-ft  depth  appears  to  be  lower  by  1  to  2 
db  than  the  noise  observed  within  1  0  ft  of  the  surface,  this  result  implies 
an  improvement  of  about  6  db  relative  to  a  single  near -surface  instrument. 
However,  this  estimate  is  probably  low  because  of  the  method  used  in 
designing  and  evaluating  the  filters.  For  the  filter  design,  signal  statistics 
were  computed  irom  entire  signal  records  rather  than  from  only  the  portions 
where  signal-to-noise  ratio  was  high.  Therefore,  the  signal  statistics  were 
contaminated  by  noise,  and  the  optimum  filters  were  required  to  pass  not 
only  signal  but  also  ambient  noise,  to  a  lesser  degree.  In  the  evaluation 
stage,  entire  signal  records  agi-’n  were  used  tc  compute  input  and  output 
nower  spectra  of  signal.  Thus,  response  functions  computed  for  signal 
were  actually  some  kind  of  average  of  the  signal  and  noise  responses.  Both 
of  these  effects  tend  to  reduce  the  estimate  of  optimum  improvement.  It  is 
possible  that  more  meaningful  results  could  be  obtained  by  employing  a 
theoretical  signal  model,  although  this  approach  has  its  pitfalls  in  that  arti¬ 
ficial  signal-to-noise  improvements  can  occur  because  of  seismometer  gain 


ine  quality. 


2,3 


The  Wiener  filters  are  designed  for  wideband  signal  preserva¬ 
tion,  whereas  the  prediction  error  filters  are  used  to  enhance  the  first  cycle 
or  two  of  the  signal.  Signal-to-noise  ratio  improvements  for  the  two  types 
of  filters  are  similar;  however,  the  prediction  error  filters  are  designed 
without  regard  to  distortion  of  the  signal. 

4 

A  deghosting  routine  was  applied  to  the  6900-ft  and  8900-ft 
records  to  separate  the  upward-  and  downward-traveling  parts  of  the  tele- 
seismic  signals.  Figure  1-12  illustrates  results  of  a?3plying  this  method  to 
the  five  noise  samples  and  two  teleseisms  recorded  at  those  depths.  When 
the  correct  time-shifts  are  applied,  a  remarkable  similarity  is  found  among 
the  estimated  upward-  and  downward-traveling  signals  and  the  signals 
observed  at  200  ft  (Figure  1-7).  Similarities  in  the  noise  records  could  be 
evidence  of  P-wave  noise. 

The  results  show  that  a  pair  of  deeply-buried  instruments  can 
be  used  effectively  to  estimate  the  signal  which  should  be  observed  at  the 
surface.  In  this  case,  the  signal-lo-ncise  improvement  relative  to  the  shal¬ 
low  instrument  is  limited,  but  the  improvement  would  probably  be  much 
better  in  high-n^isc  areas,  such  as  Grapevine,  where  the  attenuation  of  noise 
with  depth  is  more  severe. 
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Figure  1-9.  Estimated  Performance  of  the  Optimum  25 -Point  Filters 
for  an  Array  Occupying  Depths  of  200,  2900  and  4900  Ft 


Figure  I- 10.  Estimated  Performance  of  the  Optimum  25 -Point  Filters 
for  an  Array  Occupying  Depths  of  200,  6900  and  8900  Ft 
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Figure  I— 11.  Estimated  Performance  of  the  Optimum  25-Point  Filters 
for  an  Array  Occupying  Depths  of  200,  4900  and  8900  Ft 
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4.  In  vcstigations  of  a  6-Element  Arra\ 

An  attempt  has  been  made  to  evaluate  a  6-element  vertical 
array  by  means  of  the  data  recorded  at  UBO.  The  recordings  listed  in 
Tables  1-1  and  1-2  were  sufficient  to  provide  estimates  of  all  the  required 
auto-  and  crosscorrclatious  of  signal  and  noise. 

Because  of  poor  quality,  it  was  necessary  to  discard  three 
noise  samples  (recorded  on  13  October).  One  teleseism  record  (SN  18)  had 
to  be  discarded  because  of  a  computing  error  which  would  have  been  expen¬ 
sive  to  rectify.  For  each  record,  a  55 -point  deconvolution  filter  was 
designed  to  whiten  and  normalize  the  spectrum  of  the  200-ft  trace. 

This  filter  then  was  applied  to  the  three  traces.  The  purpose  of  this  step 
was  to  "equalize"  noise  or  signal  samples  obtained  at  different  times.  Phase 
responses  of  the  deconvolution  filters  were  not  determined  since  they  have 
no  effect  on  the  subsequent  results.  Aftei  deconvolution,  the  set  of  three 
autocorrelations  and  three  crosscorrelations  was  computed  for  each  record. 

Six-cl.  nnel  correlation  matrices  were  constructed  by  equat¬ 
ing  each  of  the  constituent  correlation  functions  to  the  average  of  all  the 
available  estimates.  The  number  of  estimates  varied  from  50  (1  x  1  noise 
correlation)  to  1  (2  x  3,  2x5,  2x6,  3x6,  and  4x5  signal  correlations) 
where  the  channels  are  numbered  in  an  order  of  increasing  seismometer 
depth.  The  resulting  6-channel  correlation  sets  are  illustrated  in  Figures 

th  th 

1-13  and  1-14.  The  curve  in  the  m  and  n  column  is  the  correlation  func¬ 
tion  cp  (t).  Horizontal  bars  indicate  t  ie  zero-lag  values  cp ( 0 ).  Only  half 
mn 

of  the  matrix  is  plotted  since  the  remaining  crossccrrelations  are  simply 
the  time-reverse  of  those  shown. 

For  comparison,  theoretical  correlation  sets  arc  shown  in 
Figures  1-15  through  1-18.  These  correlation  functions  were  obtained  from 

the  UBO  P-wave  and  surface-wave  models  by  inverse  transformation  of 
the  auto-  and  crosspower  spectra  after  a  scaling  operation  which  whitened 
the  spectrum  of  the  near-surface  instrument.  Agreement  is  excellent 
between  the  experimental  and  theoretical  signal  correlations  even  though 
the  experimental  curves  are  contaminated  by  noise.  Visual  interpretation 
of  the  noise  correlations  is  difficult  since  the  noise  undoubtedly  consists  of 
several  modes  with  different  spectra.  However,  there  is  a  strong  sugges¬ 
tion  that  mantle  P-wave  noise  is  an  important  consitutent  of  the  noise  field. 
This  is  especially  evident  in  the  case  of  the  1x5  and  1x6  experimental 
noise  correlations  where  spacing  of  the  peaks  agrees  with  the  experimental 
signal  and  theoretical  P-wave  correlations  but  not  with  any  of  the  theoretical 
surface-wave  correlations.  It  is  planned  to  use  a  multivariate  regression 
operation  to  attempt  to  separate  the  observed  noise  model  into  its  constitu¬ 
ent  modes. 
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figure  1-13.  Matrix  of  Experimental  Signal  Correlations 
for  6-Element  UBO  Vertical  Array 
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Figure  1-17.  Mat"  ix  of  Theoretical  First-Shear-Mode  Correla¬ 
tions  for  6-Element  UBO  Vertical  Array 
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Figure  1-18.  Matrix  of  Theoretical  Second-Shear  Mode  Corre¬ 
lations  for  6-Element  UBO  Vertical  Array 
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The  average  signal  and  noise  power  at  each  depth  is  given  by 
the  corresponding  zero-lag  value  of  the  autocorrelation.  These  values  are 
plotted  in  Figure  1-19.  Relative  to  200-ft  depth,  the  average  signal-to-noise 
ratio  improvement  which  may  be  gained  by  deep  burial  is  found  to  be  only 
about  1  db. 


There  was  an  attempt  to  design  a  6-channel  signal  enhance¬ 
ment  filter  in  the  time  domain  using  the  Wiener  -  Levinson  formulation,  but 
it  was  unsuccessful  since  the  input  correlations  failed  to  meet  the  necessary 
conditions  imposed  by  that  method.  This  result  was  not  entirely  unexpected 
since  there  is  only  a  small  likelihood  that  data  combined  in  the  manner 
described  above  will  produce  a  valid  correlation  matrix. 

Because  the  requirements  are  less  stringent  in  the  frequency 
domain,  it  was  decided  to  investigate  the  frequency-domain  stability  of  the 
experimental  model.  The  correlation  functions  were  Fourier-transformed 
to  produce  auto-  and  crosspower  spectra.  Figures  1-20  and  1-21  illustrate 
the  autopower  spectra  obtained  for  signal  and  noise  when  a  hanning  lag  win¬ 
dow  was  applied.  Frequency-domain  filters  were  designed  for  3  models: 
transform  of  correlations  with  no  smoothing,  transform  with  first-order 
hanning,  and  transform  with  fourth-order  hanning.  Even  in  the  frequency 
domain,  the  data  were  not  of  sufficiently  high  quality  to  prevent  instabili¬ 
ties.  Figure  1-22  shows  the  results  obtained,  ignoring  frequencies  where 
the  predicted  output  power  or  signal-to-noise  improvement  was  negative. 

Since  the  predicted  improvements  (where  they  could  be 
obtained)  were  less  than  those  estimated  for  3-element  arrays  (Figures  1-9, 

I  - 1  0 ,  I-  11),  it  is  concluded  that  none  of  the  results  presented  in  Figure  1-22 
can  safely  be  considered  meaningful.  The  attempt  to  combine  data  of  dif¬ 
fering  quality  and  from  different  times  has  resulted  in  a  model  which  is 
not  physically  possible,  containing  elements  with  unrealistic  properties 
such  as  coherences  greater  then  unity.  The  data,  although  they  do  not  meet 
the  requirements  for  input  to  multichannel  filter  evaluation  programs,  should 
be  sufficiently  good  to  permit  application  of  some  interpretive  techniques  to 
derive  information  about  the  generation  and  propagation  of  ambient  noise  at 
UBO. 

5.  Comparison  of  Theoretical  and  Experimental  Results 

The  design  of  a  theoretical  signal  and  noise  model  for  UBO 

has  been  described.  *  Theoretical  performance  estir  ates  for  various  array 

5 

configurations  were  presented  in  another  report.  Signal  enhancement  capa¬ 
bilities  predicted  for  vertical  arrays  from  the  theoretical  model  were  con¬ 
siderably  greater  than  those  suggested  by  the  experimental  results.  The 
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Figure  1-20.  Average  Signal  Power  Spectra  Computed  from 
Experimental  UBO  Vertical  Array  Data 
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most  probable  reason  for  this  disagreement  is  that  the  theoretical  noise 
model  is  not  an  accurate  representation  of  the  ambiect  noise  field  at  UllO. 

The  theoretical  model  was  derived  from  the  assumption  of 
uniform,  isotropic,  perfectly  clastic  sedimentary  layers  overlying  a  uni¬ 
form,  isotropic,  perfectly  clastic  hall-space.  It  was  found  then  that  over  90 
percent  of  the  noise  energy  observed  at  the  surface  should  be  propagating  in 
the  fundamental  Rayleigh  mode.  However,  if  inelastic  attenuation  and  scat¬ 
tering  mechanisms  arc  considered,  it  is  found  that  the  noise  power  become.', 
redistributed  among  the  modes.  The  data  recorded  at  UBO  support  the  con¬ 
tention  that  Rayleigh  mode  noise  is  not  the  only  significant  constituent  of  the 
ambient  noise. 


Computations  have  been  made  in  order  to  investigate  the  effect 
of  redistribution  of  energy.  Ths  signal  and  noise  power  j  icctra  previously 
assumed  for  a  near-surface  seismometer  are  reproduced  in  Figure  1-2  3. 
White  incoherent  noise  was  included  with  assumed  power  levels  of  -13,9  db 
(relative  tc  the  maximum  signal  power  density)  at  the  surface  and  -33.9  db 
at  depths  of  200  ft  or  more  (UBO  Model  20).  Two  vertical  arrays  were 
investigated,  array  B  being  identical  to  the  6-clemcnt  experimental  array 
studied  in  this  report  and  array  A  being  the  same  except  that  the  200-ft 
instrument  was  replaced  by  a  surface  instrument.  The  noise  model  has  been 
modified  by  applying  to  the  Rayleigh  mode  an  attenuation  factor  which  was  an 
exponential  function  of  frequency.  The  modified  spectrum  is  shown  by  the 
broken  line  in  Figure  1-23.  In  Figure  1-24,  the  predicted  signal-to-noise 
ratio  improvements  are  presented.  The  solid  curves  show  the  improvement 
of  vertical  array  A  relative  to  a  single  surface  seismometer,  while  the  bro¬ 
ken  curves  show  the  improvement  of  array  B  relative  to  a  seismometer  at 
200-ft  depth.  The  outputs  of  the  two  arrays  are  very  similar.  The  differ¬ 
ences  between  the  A  and  B  curves  represent  improvements  in  single  seis¬ 
mometer  performance  which,  it  is  assumed,  may  be  obtained  by  shallow 
burial . 


It  is  clear  from  Figure  1-24  that  a  redistribution  of  energy 
which  reduces  the  importance  of  the  Rayleigh  mode  is  likely  to  reduce  the 
effectiveness  of  a  vertical  array.  This  is  probably  bccaus,  the  higher  modes 
resemble  P-waves  much  more  closely  than  dues  the  Rayleigh  mode  (as  seen 
by  a  vertical  array).  Since  velocity  filtering  with  a  horizontal  array  should 
not  be  affected  much  by  a  redist  ribution  of  energy,  it  follows  that  the  pre¬ 
dicted  advantage  of  horizontal  ariays  over  vertical  arrays  should  be  increased 
when  the  noise  model  is  made  more  realistic. 

Another  probable  source  of  disagreement  between  experi¬ 
mental  and  theoretical  results  is  the  presence  of  P -wave  noise  in  the  field 
data.  Figures  1-25  and  1-26  are  contoured  plots  of  wavenumber  spectra 
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Figure  1^23.  Assumed  Input  Power  Spectra  of  Signal  and  Noise  at  a 
Surface  Seismometer,  Expressed  in  db  Relative  to 
the  Maximum  Signal  Power 


Figure  1-24.  Effect  on  Vertical  Array  Performance  of 
Attenuation  of  Rayleigh  McJe 


computed  from  data  recorded  12  October  1964  with  the  10-element  UBO 
surface  array.  Visual  examination  of  these  spectra  indicates  that  approxi¬ 
mately  1  5  to  20  percent  of  the  ambient  noise  at  1  cos  and  1.5  cps  was  asso¬ 
ciated  with  propagation  velocities  greater  than  6  km/sec.  The  symmetry  of 
this  part  of  th'  spectra  suggests  th?t  the  high-velocity  noise  is  more  likely 
to  be  due  to  mantle  P-waves  than  to  atmospheric  effects.  This  amount  of 
P-wave  noise  implies  an  ’improvement  limit  of  7  to  8  db,  even  if  the  surface- 
wave  noise  and  incoherent  noise  could  be  rejected  perfectly. 

It  is  planned  to  attempt  to  separate  the  experimentally 
observed  noise  into  its  component  parts  and  use  these  results  to  modify  the 
theoretical  model.  The  perfoimance  estimates  for  horizontal  and  vertical 
arrays  then  may  be  recomputed  for  a  more  realistic  model. 

C.  COLLECTION  OF  VERTICAL  ARRAY  DATA  AT  GRAPEVINE 

Digital  seismic  data  were  recorded  intermittently  during  tne 
period  3  March  to  16  March  1965  at  the  Trigg  research  well  near  Grapevine, 
Texas.  Instruments  which  were  employed  included  deep-well  seismometers 
at  depths  of  4500  ft  (1.37  km),  5500  ft  (1.68  km),  6500  ft  (1 . 98  km),  7500  ft 
(2.28  km),  8500  ft  (2.59  km),  and  500  ft  (2.89  km),  and  a  3-componcnt  sys¬ 
tem  at  the  surface.  The  digitization  and  recording  equipment  was  the  same 

as  that  used  at  UBO  for  vertical  array  recording*  with  high-  and  low-gain 
channels  recorded  for  each  instrument. 

Data  edited  off  the  field  tapes  for  computer  analysis  consist 
of  60  two-minute  noise  samples,  16  local  events  (such  as  quarry  blasts)  and 
12  events  tentatively  designated  as  near-regional  events  (300  series)  or 
telescisms  (400  series).  These  are  listed  in  Tables  1-3  through  1-5. 

It  is  planned  to  use  the  Grapevine  recordings  for  vertical 
array  evaluations  such  as  were  attempted  with  the  U3<”'*  data.  The  Grape¬ 
vine  data  should  be  more  suitable  for  this  purpose  since  the  entire  array 
was  in  operation  throughout  the  recording  period,  and  it  will  not  be  neces¬ 
sary  to  rely  upon  time  -  stationarity  to  construct  an  array  of  more  than  three 
elements . 
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CATALOG  OF  2 -MIN  SAMPLES  FROM  GRAPEVINE 


Date 

Code  No. 

Start  Time  (GCT) 

Date 

Code  No. 

Start  Time  (GCT) 

3/4 

500 

15:27:2  0 

3/6 

536 

22:27:10 

501 

15:39:10 

537 

22:29:10 

5  02 

15:47:3  0 

3/8 

540 

02:58:10 

503 

15:49:30 

541 

03:00:10 

504 

15:51:30 

542 

03:02:10 

3/5 

505 

15:56:10 

543 

03:13:10 

506 

16:03:50 

5  44 

03:15:10 

507 

16:05:50 

545 

03:17:10 

510 

16:17:30 

546 

1 1:45:20 

511 

21 :38:10 

547 

1 1:47:20 

512 

'  '  :00:1  0 

550 

!.  1:49:20 

513 

21  :C2 :1  0 

551 

12:07:10 

514 

21 :1 1 :1 0 

5  52 

12:09:10 

3/6 

515 

19:55:20 

3/13 

553 

00:26:10 

516 

19:56:20 

554 

00:33:10 

517 

19:59:20 

555 

00:56:50 

520 

20:15:40 

556 

05:54:50 

521 

20:17:40 

557 

06:02:10 

522 

20:34:10 

560 

06:25:10 

523 

20:36:10 

561 

06:27:10 

524 

20:49:50 

562 

06:42:10 

525 

20:51 :50 

563 

07:09:10 

526 

20:53:50 

564 

07:29:20 

527 

21:01  :30 

565 

1  0:25:1  0 

3/10 

530 

10:21 :1  0 

566 

11:25:  0 

531 

10:23:10 

567 

14:04  !  0 

532 

1  0:37:20 

570 

14:28:10 

533 

14:30:40 

3/16 

571 

1  5:54:1  0 

534 

22:17:50 

5  72 

1 6:4o:20 

535 

22:25:10 

573 

16:56:50 
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Table  1-4 


CATALOG  OF  LOCAL  EVENTS  RECORDED  AT  GRAPEVINE 


Date 

Code  No. 

Start  Time  (GCT) 

I 

Date 

Code  No. 

Start  Time  (GCT) 

3/5 

201 

20:51 :40 

3/4 

210 

16:33:20 

3/6 

202 

20:19:40 

211 

1  7:03:10 

3/13 

203 

13:49:10 

212 

17:07:30 

204 

14:41:50 

214 

17:45:00 

3/16 

205 

17:09:10 

3/5 

215 

16:55:50 

206 

1  7:12:10 

216 

1  7:01 :40 

207 

1  7:19:10 

217 

1  7:07:1  0 

220 

17:10:00 

_ 

221 

17:30:30 

Table  1-5 

CATALOG  OF  GRAPji,  .'INE  RECORDINGS  TENTATIVELY  DESIGNATED 
TELESEISMS  AND  NEAR -REGIONAL  EVENTS 


Date 

Code  No. 

Arrival  Time  (GCT) 

Source  Information  (where  available) 

Location 

Magnitude 

6  (degrees) 

3/5 

400 

21:36:23 

3/6 

401 

20:43:42 

Philippine  s 

5.8 

116 

300-305 

21:10:20 

SE  Missouri 

5.  3 

6 

3/8 

3  06 

03:08:10 

403 

12:12:23 

Alaska 

4.  5 

44 

3/10 

402 

10:26:25 

3/13 

404 

01 :3  3 :58 

Chile 

4.4 

65 

405,406 

07:42:13 

Alaska 

5.  5 

50 

410 

14:06:57 

Fiji 

5.  7 

93 

411 

14:32:36 

Rat  Is. 

4.  6 

61 

3/16 

307 

16:3  5:0- 

412 

16:59:04 

E  ol  Honshu 

5.6 

8  9 
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SECTION  II 

DECONVOLUTION  AND  AUTOCORRELATION  STUDY 


The  application  of  deconvolution  and  autocorrelation  techni¬ 
ques  to  three  ensembles  of  earthquake  records  from  the  multichannel 
processor  at  CPO  is  being  studied.  Statistics  of  discrepancies  between 
observed  and  USC&GS  reported  origin  times,  magnitudes  and  depths  of 
focus  are  examined. 

Ensemble  and  subensemble  average  autocorrelations  are 
presented  for  Ensemble  I  (73  Kurile  Islands  events).  Very  little  difference 
can  be  observed  among  the  subensemhles  or  between  unweighted  and 
normalized  averages.  It  is  hypothesized  that  certain  features  in  the  auto¬ 
correlations  are  due  to  reflections  from  the  Conrad  and  Mohorovicic 
discontinuities  under  CPO. 

Ensembles  II  and  III  (60  and  119  events  of  worldwide 
distribution)  are  described  and  the  results  of  individual -event  decon¬ 
volution  presented.  Although  signal  waveforms  are  contracted  somewhat, 
there  is  no  obvious  improvement  in  detectability  of  later  phases.  The 
nature  of  the  ambient  noise  recorded  at  CPO  is  discussed  briefly. 

A.  DESCRIPTION  OF  DATA 
1.  Introduction 

Three  distinct  groups  (ensembles)  of  earthquakes  recorded 
at  the  Cumberland  Plateau  Seismological  Observatory  (CPO)  are  being 
studied  by  use  of  deconvolution,  autocorrelation  and  depth -of -focus 
estimation.  Analog  FM  tape  recordings  of  the  output  of  the  CPO  array 
processor  were  converted  to  digital  format  for  computer  analysis.  En¬ 
semble  I,  consisting  of  a  suite  of  73  events  from  the  Kurile  Islands,  was 

de  scribed  in  detail  in  an  earlier  report.  Ensemble  II  consists  of  ail  non- 
Kurile  Islands  USC&GS  reported  events  between  25  September  1963  and 
18  October  1963.  Ensemble  III  consists  of  randomly  selected  events 
reported  by  the  USC&GS  between  20  February  1963  and  14  May  1963. 

Figures  II - 1  and  II -2  (CPO  polar  projection  maps)  show 
the  random  geographic  distribution  of  the  86  Ensemble  II  events.  Figures 
II-3  and  II-4  are  similar  presentations  for  the  174  Ensemble  III  events. 

The  multichannel  filter  operating  at  CPO  when  Ensemble  III 

7 

was  recorded  was  IP -8.  Ensembles  I  and  II  were  recorded  while  pro- 

,  .  8 
cessor  MIP-4  was  m  use. 
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Figure  II- 


3  r» 


2.  Geographic  Distribution  of  Ensemble  II  Events  Beyond  90° 
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Figure  II- 3 *  Geographic  Distribution  of  Ensemble  III  Events  Within  90° 


II-4 


2.  Arrival  Times 


Observed  P  arrival  times  for  Ensemble  II  averaged  0.  1 
f  1.9  sec  later  than  those  predicted  from  Jeffreys -Buller.  (J-B)  travel¬ 
time  tables  and  USC&tGS  reported  epicenters  and  origin  times.  Ensemble 
III  observed  P's  averaged  0.  !  ±  1.7  sec  ahead  of  prediction.  Figure 
II-5  is  a  normalized  distribution  plot  of  observed  -minus -predicted  arrival 
times  for  all  three  ensembles. 

For  some  distant  events,  PKP  was  the  first  arrival.  In 
Ensemble  II,  observed  PKP  averaged  1.  3  ±  1.9  sec  later  than  predicted, 
while  observed  PKP  in  Ensemble  III  was  0.  3  +  2.  6  sec  later  than  predicted. 

The  curve  for  Ensemble  I  in  Figure  II - 5  shows  a  bias  of 
about  +  0.  4  sec,  whereas  Ensembles  II  and  III  show  a  much  more  even 
distribution  of  arrival -time  anomalies.  This  bias  could  be  due  to  some 
properties  of  the  travel  path  common  to  all  events  from  the  Kurile  Islands. 

q 

These  results  differ  from  those  obtained  in  another  study  where  average 
time  residuals  for  CPO  were  determined  to  be  0.02  +  1.0C  sec  for  270 
Kurile  Islands  events  and  -0.  26  ±  0.  98  sec  for  609  events  with  worldwide 
distribution. 

3.  Magnitudes 

Magnitudes  computed  from  trace  measurements  using 
prescribed  USC&tGS  techniques  for  Ensemble  II  averaged  0.  3  +  0.  04 
less  than  reported  by  USC&iGS.  Those  computed  for  Ensemble  III  were 
0.  1  ±  0.  5  less  than  reported.  Figure  II -6  is  a  normalized  distribution 
plot  of  observed-minus -published  magnitude  for  all  three  ensembles.  The 
curves  for  Ensembles  I  and  II  show  more  of  a  trend  toward  negative  differ¬ 
ences  than  Ensemble  III.  One  reason  for  this  difference  is  that  a  narrow - 
band  filter  was  applied  to  Ensembles  I  and  II.  This  filter  had  the  effect 
of  lowering  the  trace  amplitude  and,  hence,  the  observed  magnitude. 

4.  Depth 

Observed  depths  of  foci  were  calculated  from  pP -P  arrival  - 
time  differences.  Figure  II -7  is  a  normalized  distribution  plot  of  observed 
depths  minus  published  depths  for  the  three  ensembles.  As  seen  from  the 
figure,  100  percent  of  the  events  for  which  pP  was  observed  in  Ensemble  I 
lay  within  th  '  +25  km  accuracy  professed  by  the  USG&cGS.  Ensemble  II  pP 
depths  showed  91  percent  agreement  with  the  puolished  depths  +25  km 
while,  in  Ensemble  III,  88  percent  of  the  depths  determined  from  pP  lay 
within  +25  km  of  published  estimates. 

Figure  11-8  is  a  histogram  of  pP  observance  over  the  magni  - 
tddo  and  epicentral  distance  ranges  encountered  in  each  of  the  ensembles. 
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Figure  II-7.  Cumulative  Distribution  of  Differences  Between  Depths  of 
Foci  Computed  from  pP  Arrival  Time  and  Estimated  by 
USC&GS 


It  is  worthy  of  note  that  all  of  Ensemble  I  (the  Kurile  Islands  data)  would 
necessarily  occ:r  about  85°.  There  is  some  evidence  from  this  figure 
that  intermediate  magnitude?,  in  particular  5.0  to  6.0  magnitude,  yield 
the  highest  percentage  of  pF.  The  pP  phases  from  large  events  often  are 
masked  by  the  P  wave  train  and  pP  phases  from  small  events  often  are 
hidden  by  noise.  The  optimum  magnitude  for  pP  observance  is  generally 
less  in  the  case  of  near  events  than  more  distant  ones  because  the  P  wave 
train  from  near  events  is  more  elongated. 

B.  AVERAGE  AUTOCORRELATIONS  FROM  ENSEMBLE  I 

In  an  earlier  report,  ^  use  of  deconvolution  filters  designed 
from  ensemble -average  statistics  was  described.  Results  were  found  to 
differ  only  slightly  from  those  obtained  when  each  record  was  deconvolved 
using  a  filter  based  only  upon  its  own  statistics.  In  that  study,  the  ensemble  - 
average  statistics  were  determined  by  computing  an  unweighted  mean  of  the 
autocorrelations  of  the  7  3  events.  It  was  suggested  that  a  method  should 
be  investigated  in  which  the  events  are  normalized  to  the  same  mean -square 
vaiue  before  averaging. 

There  are  advantages  and  disadvantages  to  both  approaches. 

If  an  unweighted  mean  of  the  autocorrelations  is  used,  only  the  large  events 
will  have  any  significant  influence  on  the  results.  If  the  autocorrelations 
are  weighted  to  normalize  their  zero -lag  values,  added  importance  will  be 
given  to  small  events  which  generally  have  small  signal -to -noise  ratios. 

Thus,  the  output  from  the  normalized  set  may  be  expected  to  be  contaminated 
by  noise  to  a  greater  degree  than  that  from  the  unnormalized  set. 

Unnormalized  and  normalized  average  autocorrelations 
computed  f^  r  six  subgroups  chosen  from  Ensemble  I  and  for  the  complete 
ensemble  are  illustrated  in  Figures  II - 9  and  11-10.  Autocorrelations  were 
computed  to  lags  of  24.8  sec  using  45-sec  records  starting  just  prior  to 
the  P  onset.  The  groups  represented  in  Figures  II -9  and  II  - 1 0  are  as  follows 

(1)  31  events  with  observable  pP  and  a  uniform  distribution 
of  P -pP  intervals 

(2)  15  events  with  large  amplitudes 

(3)  26  events  with  medium  amplitudes 

(4)  30  events  with  small  amplitudes 

(5)  30  events  with  clear  later  phases 

(6)  42  events  with  "ringy’’  autocorrelations 
(this  phenomenon  is  apparently  an  effect 
of  source  mechanism  since  it  does  not 
appear  consistently) 

(7)  All  73  events  of  Ensemble  I 
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Figure  II-9.  Unnormalized  Average  Autocorrelations  of  Subgroups  of 
Events  from  Ensemble  I 
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Figure  II- 10.  Normalized  Average  Autocorrelations  of  Subgroups  of 
Events  from  Ensemble  I 
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It  is  difficult  to  find  any  meaningful  difference  either 
between  the  normalized  and  unnormalized  autocorrelations  or  among  the 
autocorrelations  of  the  various  subgroups.  For  this  reason,  it  was  con¬ 
cluded  that  probably  no  advantage  would  be  gained  by  computing  additional 
deconvolution  operators. 

There  is  clearly  not  enough  resolution  in  the  autocorrelations 
to  permit  conclusions  about  near -surface  layering  at  CPO  There  is  a  sug¬ 
gestion  of  persistent  features  in  the  autocorrelations  at  times  of  approxi¬ 
mately  6,  10  and  13.  3  sec.  These  may  be  associated  with  primary  and 
multiple  reflections  from  the  Conrad  and  Mohorovicic  discontinuities  under 
the  receiving  station.  Additional  statistical  significance  will  be  provided 
by  the  results  for  Ensembles  II  and  III  to  test  this  hypothesis. 

C.  INDIVIDUAL-EVENT  DECONVOLUTION  OF  ENSEMBLES  II  AND  III 

From  Ensembles  II  and  III,  a  selection  was  made  of  60  and 
119  events  with  record  quality  sufficiently  high  to  permit  useful  analysis. 
Figures  II - 1 1  and  11-12  contain  playbacks  of  the  selected  events  with  no 
differential  magnifications  applied.  In  order  that  the  signal  waveforms  may 
be  easier  to  observe,  Figures  11-13  and  11-14  present  the  same  events  with 
magnifications  applied  to  make  all  P -wave  amplitudes  approximately  the 
same.  For  each  event,  a  55-point  (5.  5-sec)  dec  nvolution  filter  v/as 
designed  to  whiten  the  spectrum.  The  results  oi  filtering  the  records  with 
their  associated  deconvolution  operators  are  shown  in  Figures  11-15  and 
11-16. 

Source  information  for  each  event  is  presented  in  Figures  II - 1 1 
and  11-12.  To  permit  a  longer  interval  to  be  shown,  this  information  is 
omitted  from  Figures  11-13  through  II- 16  which  display  the  events  in  the  same 
sequence  as  Figures  II -11  and  11-12. 

Playbacks  cf  the  raw  data  (Figures  II - 1 1  and  11-12)  provide 
valuable  information  regarding  the  nature  of  the  ambient  noise  at  CPO  and 
the  relative  merits  of  the  two  multichannel  filter  systems.  Ensemble  III 
contains  recurds  output  by  processor  IP-8,  a  system  designed  from  a 
the  or  ctic  al  model  of  ambient  noise  assumed  to  be  in  the  speed  range  of 
2.  5  to  3.  5  km/ sec.  Subsequently,  the  theoretical  model  was  found  to  be 
seriously  in  error;  noise  velocities  as  high  as  4.  5  km/sec  were  measured. 
Multichannel  filter  systems,  such  as  MIP -4  which  produced  the  Ensemble  II 
data,  were  designed  from  measured  noise  statistics.  Also,  low-cut  filters 
were  added  to  attack  the  6-sec  microseism  peak.  Since  no  differential 
magnifications  were  applied  in  preparing  Figures  II - 1  1  and  11-12,  the  only 
trace -to -trace  inconsistencies  in  amplitude  should  be  due  to  adjustments 
made  to  the  recording  and  processing  systera  at  CPO,  and  it  is  understood 
that  these  effects  are  of  a  minor  nature. 
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LKs^ _ 


(i2a) 


(12b) 


3-1 

3-2 

3-3 

3-4 

3-3 

3-6 

S-Y 

3-8 

3-9 

3-10 

3-11 

3-12 


Event 

No. 


■  (  « ) 


204.2 
19?.  1 

200.  6 

302.0 

301.0 
134.  y 

19?.  8 
140.  3 

190.4 

202.  1 

196.0 

220.9 


Focal 

Depth 


130 

150 

33 

33 

>3 

33 


eiiliiS 

'yr'h* 

vVy*V/vs^v\  /CvV^'  V* 

*  •  i  *  it  # 


33 

33 

33 


%.  1 

4.6 

5.5 


!*  ,V»Vi 

.  '  ,«  *  *1  |  4  *  < 


■»  ■  ,*'Y 


•  s  '  k  ,  A  .  .  *,  It  ' 

\  \!  •  \\  >  f  J  ■  ■  .* 


..  i  - - 

4. 


3-13 
3-1  4 
?  15 
3-16 
3-17 
3-18 

3-11 

3-20 

3-21 

3-22 

3-23 

3-24 


M  V» 


22.63 

22.88 

23.0? 

23.  15 

23.  15 
23.41 

23.76 

23  87 

25.94 

*6.02 

27.  15 

28.53 


193.  8 
219.  3 
129.  8 

18*.  5 

184.5 

186.4 

279.6 

185.4 

**.  5 

239.6 

278.  2 

172.9 


33 

33 

39 

34 

33 


44 

16 

53 


«'l  *  ' 

3 '  7  s>m*t  — 

<•■  v-w^,  ':$fJp:rpfy ,  Vv 


4.  b 

4.3 

4.0 


■w-X/.VV  /  \ 

•-‘  i  f  i  n  i  mi  mmmfflTf  mnmniiiTTtiTiTtitrmnTi 


Timing  Lines  -  1.0  Sec 


(I2d) 


(12e) 


Event 

No. 


1-17 

1-18 

1-1? 

3-40 

1-41 

1-42 

1-41 

1-44 

1-45 

1-46 

1-47 

1-48 


Event 

No. 


1-49 

3-50 

3-51 

3-52 

3-51 

3-54 

3-55 

1-56 

1-57 

1-58 

1-59 

3-60 


36.95 

17.32 

17.61 

39.3 

39.4 
40.06 

40.06 

42.08 

46.  C6 

46.44 

46.52 

47.26 

t 


168.6 
16  7.  1 

100.  5 

107.8 
108.2 
173.  8 

171.6 

173.1 

32.36 

32.35 

170.4 

12.65 


Focal 

Depth 


47.28 

48.  14 

48.  18 


49.17 
49.5 

51.17 

51.  35 
51.92 

52.  15 

51.4 
54.  1 


126.8 

29.5 


89 

105 

13 

33 

33 

31 

31 

33 

61 

78 

13 

39 


Focal 

Depth 


4.  9 

4.  5 

4.5 

5.4 


*  <  ,  »*  *  *  * 


.  VMvMWW 

i  V  il 


29.5  *  15 


29.0 

215.4 
116.  9 

23.0 
164.  7 
161.8 

163.  7 

8.7 

116.0 


5.  7 

5.0 

5.5 

Mag. 


(  i 

mtTtnmfj  ifi 


no 

105 


33 

31 


4.5 

4.6 

6.  0 

4.  6 

4  6 

5.  1 

4.6 


.  1|  t  ,  , 


•  ^  i  ’  1  ,i  t 

iffy* t 

A,*  ’ / *, 


54  ^/''y , 

5- 3 


4.6 


Figure  11-12. 


IP-8  Output  Records  of  Ensemble  III  Teleseisms 
(II-12a  through  12e) 

II-  17/18 


(12f) 


(12g) 


Ev«nt  .  Foc^  w  rrr 

No  &(*«)  Aaimuth  Depth  Mag.  j 


Event 

No. 


3-73  63. 


\  y  k  J f  ^ 

«  «•«  /\JXT\/\  C  AA/V 

202  6.0  ' V--*— ^  |  f'/  \  \ 

13  5-0  VVXVV1,  '  ^  '•  ’  '  -  '/  '  7/  /  A‘ 


U  6.8 

271  4.2  W"  \<V\V 


~.  J>  /V,VS/VVV->VV*^^  V  A’,', 

J  i  1 

-V\\  v1'(*^vw^‘,',V,V,  v  y^W^n^.w/v^W.’.V.vww 


Focal  | 
Depth  |  M* 


15  *  316.4  33 


3-74  6b.  15  160.  5  32 


.  53  316.  9  50 


s-)\/V^  '\'  •  ,/  /'  ’V;  '  '.A" 

4. 9  «***« —  A4wyA^*/vyv  yWstjftsY*' 


.72  3 lt>.  5  33  4.1  - - 

.62  206.8  »  5.6 

.71  206.2  33  4.9 


.62  206.8  33  5.6 


3-79  68.85  163.9  232 


3-80  69.  14  164.4  33 


.31  164.5  33 


.  38  171.  8  33 


.83  326.4  33  5.  1 


3-84  72 


194.2  33 


tmtmttmnmmffmim 


6  (***)  I  Aelmuth  j  Depth  Mag. 


(12h) 


3-86  74.5 


3-87  74. 95 


3-88  76. 13 


3-90  78.07 


3-91  78.24 


3-92  79.51 


3-93  80.48 


3-94  81.16 


!l!llll!llUliillM!iiMW 


I  IINUl 


13  4.7  j- yA'v W**'  v^/wv^/Vs/S/^ 


.  4  ✓V’W'w  WS/A 

.  4  vAa-A  ■m*’ 0«w’' \V  ^ 


i  '  '  ,  /  * 


Timing  Lines  =  1.0  Sec 


(12i) 


(12j) 


Figure  II-  12. 


Event 

No 

£('■•) 

Aiimuth 

F  oca! 
Depth 

Mag. 

3-97 

83.0 

323.0 

33 

5.2 

3-98 

83.  3 

54.  5 

33 

5.0 

3-99 

83.  31 

54.  2 

33 

4.5 

3-100 

83.43 

54.3 

33 

4.4 

3-101 

84.63 

50.  1 

42 

4.  8 

3-102 

86.46 

3267 

100 

4.  7 

3-103 

87.  33 

325.4 

110 

5.6 

3-104 

88.49 

329.6 

31 

4.  8 

3-105 

89.  5 

328.8 

255 

5.6 

3-106 

97. 9Z 

37.2 

33 

5.2 

3-107 

99.07 

327.2 

33 

5.  9 

3-108 

108.25 

246.2 

45 

7.0 

Event 

No 

Ma,»> 

I 

Aiimuth 

Focal 

Depth 

Mag. 

3-109 

108.  56 

246.0 

50 

6.2 

3-110 

120. 38 

331.2 

51 

5.0 

3-111 

12o.6 

320.3 

51 

5.  2 

3-112 

126.87 

319.0 

117 

5.0 

3-113 

127.6 

298.  8 

40 

5.5 

3-114 

130.42 

300.  7 

31 

5.6 

3-115 

137.61 

222.9 

39 

5.0 

3-116 

141.  71 

305  8 

33 

5.2 

3-117 

143. 84 

356.4 

30 

5.6 

3-118 

145.  72 

310.  1 

33 

5.6 

3-119 

148. ?4 

335.9 

33 

5.2 

s  7  V/v'M 


/  ‘  iwwMfWWrfi 

- ‘i .  W{\  Z/',* ' V  7W\‘A 


//, "AW-v\\^nvvv. 


,  yv\xvV/\A<'«W/^.'' Aa*/v'*-'**>vWW^Wv 

iifiminniTitinimrrtitffltittinmitiTmTfmri 

ifl 


,,<•.'  I  1 


^y\>  ’ ;  ( * 

*’  j 1  yW  'v ;  V  7  ■  ;  V  W* 

*  *  •»♦'  i*  •  *'»  * 

S.6  .****■<— ^«V\A<~^»vV»V'  WWr  V  wVSjfiSW** — Vv/W^W'-*^** 

*  I 

mtTmrtttmfiTntTTmTrTTmTmnmmmmnn 


IP-8  Output  Reco  'ds  of  Ensemble  HI  Teleseisms 
(II- 1  Zf  through  12j) 


TT-1Q  /  20 


(13a) 


(13b) 


(13c) 


Evont 
No. 


2- 

2-2 

2-3 

2-4 

2-3 

2-6 


I  I  I 

j  I  • 


2-10 

2-11 

2-12 

Evont 
No.  | 

2-13  i 

2*14  . 

2-15  . 

2-16  „ 

’-17  . 

2-18 


— !*' yJt  /f ;  /.WY, /V/V,  V 

*i  •  ,  i  ,  , 

-a,  yi  V,W,’,'Wiw»V,  V\,\  ( ,^V\WI Vi  (*'i'wA‘  jA/a/V'~wws«vv 

\  \^Vj'v^.Mi^i/--v^iuivyw,/~\'^-^\vwwAyw~/^AyvN/vv\^-^ 

* 


III  1  1 

111 

J  fill 

,—11 

Wijp  i  1 

*1  ,  , 

-vV  V'i  •  "VS,  .'S’;  Vvv>w/,> 

I,  •  <  *  * 


y /v/.W'  ■ 

—— — '  ^V*  I 

2-24 


imuiimiium 


2-2; 


ini 

I  i  i  1 

4^W\ 

I 

iil! 

ill! 

mr 

hi 

i."  «v  ^A/S>. 

i  mnur 

PH!  !i  llllllllil 

2-26  v^/V-vV^ ./VSf’W'W  V  •// V  (%'-v^AT  ,‘1'.’^’/.  AV  v.v/ * 

2-28  — ■ — .WW'.y^A  ,' V,'/vV,W,N 


2  29 
2-30 

2-31 

2-32 


t ‘^VwViAAAmv'  ' V|/V);V* * \fr^ l /*}+}& 

-j  '/gvW^/V^V-i^VA  //,v/V  ,V,  V. .'  j^V-AA  /.V^a'a 

ji  '.^ /  ‘ '  1  V«  /A'(V'i'W«/.'JAW,'''/-'',,,,Jv'*,'v'^vv 


2-  34  i^/^\*/i?vvvAX/(/'^W' <V\V' ^W\^^VvN^AV^vV'^fV--V/^V\\^^Av*N,VvVWVV- 
-i  V1v^1*A^\v/v.^WvA^<sV/vwyvi'v 


Timing  Lines  =  1.  0  Sec 


(13d) 


(13e) 


Figure  II-  13. 


Event 

No 

2-37 

2-38 

2-3  <* 

2-40 


I  I  | ;  i '  i I •  r . : I :  |  I  ,  I  '  j  j iFliTi  I ! 

u  ;  1 1 1  j ;  1  ii  i  1 1 1 : 1  ■ '  Mil 


-Av,  •'/tSyssi  ‘  ,V‘S'.  V 

2-41  vAv-^»vvv^^^WA  *  v ■■■  . . 

\  v  /wV'  Vs^M’.V'wvvrf/  nVV/  V-  V,  AV  (V\A7v*"v'-tyw'/ W\A~VWVr 


2-43 

2-44 

2-43 

2 -4b 


f'/'.Vl'.Vj'V  (*,vW'bY,‘Ay^NM/wCVAA/*^~^Jw>vVs~ 


’a.Vv.'^waV^SV  Vv-v'.vM  vvva»,  7 


- — -A*  ’  ^Va  a  ,V 


•“II— 


Event 

No. 


IFlIiiijii^ 


•  *  ? 

l,r-'^w^vs /-s/v*  j* 

\a/s#vwva/v*vv^V\ 


2-51 


2-52 


2-53  '^*n/WV\A^  //Vv^s^ 

i  « 

2-54  vvYv^^AWs'(  Vk 

\  ,  .VAaVi 

’  ,,yv';/A’v^^vv\r  V*;  /7/^V 


2-55 

2-58 


2-58 

2-59 

2-60 


VVV*  /.Vaw\^ 

>vw\~-i-vW;’  »  ‘,V,'V,  rfVtV ,Ws^\VWAtVW>\\V>-’ 


rn*  1 1  ftttTtT^ffifirtTt'nttnii'iTTT^^ 


Ensemble  II  Teleseism  Records  with  Differential 
Magnifications  Applied  (II -l  3a  through  12e) 


TT-?  1  /?? 


(14a) 


(14b) 


(14c) 


Event 

No 

3-1 

3-2 

3-3 


■  ' )'  ‘/.'Vi  W;.'; 


V.vV/VsA^-^-vVjV/^iVA,— AV>AAi>"A/»**W1V 

^Ay/wV/w^V ftl:*'/-*  '//■  yv.’A*',"/,* 

;*.fvs '/'//  Wv, 

-*•**#}*&**/*>  ‘ V ,‘ *■' * '  /Wv.  V.A'V/.'V/ v\ v/v^ 

\\‘i‘  vs,’  '^J'xjy' 


^vwyv**/w,V  AA.VV  >'•  v,  •,V/1,,V.VIV^*1  ,v^/VA',vvvfA W,niV#/VvVWv^' 

5-12  rtnntirrtrrnitt^  nw-vrir^^mt^ 


Evidi 

No. 


3-13 


[II 

n 

ir 

IT 

rn 

un 

iTTTTTTrTn^nirrnTiiTi 

ii  liiilii  'U 

iTTirniii!  FUTTiTn 

jj_ 

jjl 

II  1 1 1 1 1 1 1  i  1 ! 

3-14 

3-15  v^/ 
3-14 
3-17 
3-18 


*, 


•  v,‘V^v‘vV 


■  i^Mi  V,wj'v^v'v*v\  ^-^v- 

* ^\yv* /W^W/^//W -WW vV/^WV^r^/v^^^v,W^Y 
«  * 

— ;  ^  w^wvn^w 

—4  V'JWv-v^w-vViS  \*<V  yv  *,■%  aA /-v-^^vvkaiv-v- 
*  <  * 

Wv^V  A  ■'M a\W/aVvsVa^w^^ 

'Vs-/^ 


3-19 

3-20 

3-21 


3-23 


*•■<  tttrcTmttrmttttttf^^  ifflfflnti 

nnriiiiiiii! 


^  **;  y»\V 

v> 

*1  * 

3.33  4,  11.  ■  ■■ 

J-34  >A>«\^wVvwk‘.,(Vv^//  A/Wv>-vS^ 

■•“  mritmim^ 


Timing  Lines  =  1.  0  Sec 


lliillll 


uirrnrr 


'iurnn  !ii  I 


III!  M 


(14d) 


3.3*  W'^VA-A''''A/ 

3.41  *"*■■***-,< — »A^'<**— *• 

3-42  - - - <  '*  \  /V  y\OWW'^^yVVs/V'/\^''^^ 

3-43  - - - - - - - *  ^'‘JV/1V/,'V  /^/WSN^/WVn/VWS^V*  -VSA~- VvV4 

3.44  -  . -  .  - . . . —  ..  -  ■  vV^^MA-V'y"/V^AA/\A/\^wS/w'<^/w^'~ - — 

f 

3-46  !  !fa‘etor*Wf\j  ‘vWWAV'f 

>-“  IT!  I ,  tmtmtntnt^ 


iiiMiifiiiinwiiniiiii 


*A»V\  ^  VWiWl/W^V^'W* « »■ 

"'’■  '  ’  v  r  .  •  , 

v'V  ’  Y» h\  ' 

*  *  1  ,  i  r 

*<//  j  ’.Vs 


U4e) 


Ys  /'/W  ’//vV  Y/AV^V^r-VV^v//\A^V^0  vV'-'WVWW- 

*  ‘  ’•  .  . 

fj+sytiJ^  /-  v/w/^ww /vvy/'vM^^^-vv^o.~vA>,,/w^A'~~\~' 


y,^v'J*.\’'rJ*>V, 


/v^J*vv^'1V,ilyvV( *v*  'y  a\^w\VaV^a>'VV^j^VS/vVW 


f.  ft 


Figure  11-14.  Ensembie  III  Teleseism  Records  with  Differential 
Magnifications  Applied  (II-14a  through  II- 14e) 


11-23/24 


<14f) 


(I4g) 


(14h) 


!iir  nrmninT i 

iir  ill]  niiimi 


3-67 

3-68 

3-69 

3-70 

3-71 


My/-<w — -w - ^v-v^*  ^/,l'y^-WvVvA^'V^'^vA'' — Z*' 

-/  •’  V  'W V  Ita  /\  v  / V/,^/^W^>vvV 
'*  ,  ’  »  '  »  «’  »’  1 

-V/  J^Al',,Vw'V^^■^/V^V''■A/'v 


^VfwAy. 

i  (  , 

— —  ^ ' yf*wvvv*++A/v*\+*^^  /^aT^***  V 

— — \,(N^VA.V,V,V#  V  V^/^V-VW^VA/N^^VA^ 

I*  i  *  ,  * 


3-72 


■  mttTmtrnmftfitt^^ 


:„  MMMMimiiniiiiMniiji 


Ullll 


’  i/Vy^  AV— 

'^ww^'v',vs,'(v/ v  *'/V  ,*,/  V  •ys',  »V.V‘*  •’•’»'*  V 

«'v/vv7i v> 


3-83 

3-84 


/^<w, 

1  \  vNav^VvV^v/ /'vvvvZ/V/' 'v.*1*^'/**'  ’■'** — *\wuv- *^vvvv 


... niiiiiijiiiinnijiijiylimii  lllilili  iliilittliflillilf 


•  y/  *S\jW /%’-^~^v/»/'V^vA7^JW^r ' W'v'M'W' — Vj — \a/^/\J\ 


zZ/V^Aa^Aa/Vv'Vvva/^^  /-v 

•wvsAjvZ  ^  '  [ff  V  W/AW/V.'V/  yV'A'— V-'  ^jV>- 


3-94 - - 


Timing  Lines  =  1.  0  Sec 


(14i) 


(14  j) 


Figure  II- 14. 


3-103—- 

3-104- 

3-105  - 

3- 106  * 
3-107' 


i  S  I  j  !  I  i  i  !  Li  i  I  I  i 


IIPM0WPI1RW! 


',  V  V/\f  \  ;,\  r' j  ,^iYi'c/sA\Y/^V'^V'^v'<f 


//WV*A  V/>\^wAAhvwv\/^*/\A^«/^V^WV \^WV^  ^ 

y  *  1 

y,V/'V\v^wy,^,WvW 

*  .  , 

'Vvv yyv«^wy\ !^\n wA^awwva^^' vvWv/^a^n/w^*^ 

>-  ■“tTTrmTfrrtimrmiTT?nTrrm  iHmh  h  ;  n  i  iti  i  ft  n 


3-109 
3-110 
3>  111 

3-112 
.3-  a  1  3 
3-114 

3— 115  ******* 
3-116* 


i![ij!iiiiiii[iiiiinniiiiiifiiiiiiffl 


3-117  . 


/V/^kv> 

,  »  >  ,  . 
;v  /j^'Wvvy/V .  y.^\v  vWVWi/Wty 

aV*/M^  y  f,  "  fyw 

I  I 


—— ^  V  y  ■  «■*  ■!  %!.m  i  ,  i 

i  *  *  * 

— — Y / ///VW'<,V  y, 

'-w\AA-^-vVyv  WAV/  ^ 


3-119' 


3-120 


IWItniTM 


Ensemble  III  Teleseism  Records  with  Differential 
Magnifications  Applied  (II-14f  through  II-  1 4 j ) 


H-25/26 


(15a) 


(15b) 


(15c) 


2-4 


/  VvV;' 

/,  vAr^wv^~  V  \w// y^\«\*v 

t 

2-s  /n/-s^v'V^/«\A-'v\ 

z-6  \~-v- — — - ^  /^\NMn/V/vW^('^ Lf^S  .'A^W 

2  -  7  ^\/ww-v — w-^vvvv  ^v/4v,  VVAir^ V,^  ^v\y//V,V^,  yyA — 

2-8  ~A*vv/\/,*,*'VW,v/\^y^l^^  y4\  /^tyv**v,V\r'‘VsVV'Av«*AvvS/NAr^ 

2-9  Vvvwvv  Avs/*-v/|  ly*A*YAAtyVWWAy''^.^vA'VV^ 

2  - 1  - — - i\y,v,',V  V.'i,/A^Vv,VA^A,.WjV/i</,^VX^V^^V\(VtAAA'^ 

2_  1 1  yVV*^*^\A/V-'^'^wy%/'|l1Vv  *^X/Vyl»\V,V>*^V'Vi^^^■yw,■,*“''”'yA/''~',/VV^fwVVW-^ 

>••■  ittttmmTrmmtvH^^ 


Event 

No. 

2-13 

2-14 

2*15 


■ 


;;s^wwv 


V#WW»w\ 


I 


— V  //  ’v \V  y,’  *»#  \i\\l*y\  .Vv^/,  / ,WV/WW^*^ 

•  *  I  *  *  1  ' 


-  •v'f’fjf lyy'y»^Awwv;V'A^r-V-v 


2-»6 


2-17 


»  *i 

2-17  - - - - /Srf** 

2-22  -_V6 /V\A-/-^/vsyv'^VWlV*AV^VAv','^v/,*'^ 

rnnwifm^^ 


2-«;4 


Event 

No. 


2-25 
2 


I  ! 


2-30 


26 

2-27  v' /vV-^'AA^a.;^~' 

2-29  ^VvV\/'^-^^Ayvvv^1  /vW/‘yvW/'V/v«V/Iv^vv'V^/M f^/vvyvA^r 

■y^vw 

yvN^V  '//Lf  W  •  yty/w*' 

2  -  32  .^vv^v^AM/vvV-yyy*^  \\''WWV*V^ty*jAfS^^ 

2-33  ZvV^v\//^v//\AV'/^-/-  ,V/ A/^^A/J^/v^VS^V^v^'W^"v''vv,'’^VA/AA*/kVV^A/',J',U0“W^v"W 

.  V^Vyw^^  WW>'-'^'/<V*^W'‘  ■‘w‘^^-^sy'A^^^^*W^yVV*\A-Vv2/' 


2-35 


2-36 


Timing  Lines  =  1.0  Sec 


(15d) 


(15e) 


Event 
No 


2-  39  — — - — ^-vvv^. - v*v  A >_»A  ‘.‘, ;  *//.VV-  v"/A'vn\''JrvV/ ^V^V^W'vV^^VV* 

2.41  \  - — • — v^-AA/^W^/  V ‘ — . . .  . .  — — 
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Figure  11-15.  Ensemble  II  Individually  Deconvolved  Events 
(H-15a  through  15e) 
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11-16.  Ensemble  III  Individually  Deconvolved  Events 
(H-l6a  through  l6e) 
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Figure  II- 16.  Ensemble  III  Individually  Deconvolved  Events 
(II  •  I6f  through  II-l6j) 
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It  is  clear  that  processor  IP-8  permitted  a  large  fraction  of 
the  surface-wave  noise  to  pass  into  the  output.  The  level  of  the  noise  is 
found  to  vary  over  a  range  of  more  than  10  db  from  trace  to  trace.  The 
high  noise  power  observed  on  some  records  suggests  microseisrr ic  storms 
of  intermittent  occurrence.  By  contrast,  the  noise  observed  on  the  Ensemble 
II  records  is  remarkably  time-stationary,  it  is  hypothesized  that,  because 
MIP-4  and  the  low-cut  filters  were  so  much  more  efficient  at  rejecting 
surface-wave  noise,  the  noise  which  was  passed  was  composed  mainly  of 
P-waves  and  the  intensity  of  this  background  of  P-wave  noise  did  not 
fluctuate  much.  These  results  imply  that  a  detector  employing  decision¬ 
making  equipment  should  be  able  to  operate  effectively'  on  the  output  of  an 
array  processor  such  as  MIP-4  which  is  capable  of  rejecting  surface-wave 
noise  efficiently. 


In  Figures  11-15  and  11-16,  it  can  be  seen  that  individual- 
event  deconvolution  has  succeeded  somewhat  in  contracting  the  waveforms 
of  the  P 'events.  In  many  cases,  the  direction  of  first  motion  is  easier  to 
determine  since  more  energy  has  been  concenti  .ted  in  the  first  half-cycle. 
However,  many  signal  waveforms  remain  extremely  complex.  It  would  be 
expected  that  deconvolution  should  also  contract  the  waveforms  of  the 
later  arrivals,  but  these  are  no  more  apparent  on  the  deconvolved  records 
than  the  raw  data.  Presumably,  this  failure  is  caused  by  differences 
between  the  waveforms  of  P  and  the  later-arriving  phases. 

It  is  planned  to  compute  autocorrelations  for  each  of  the 
events  in  Ensembles  II  and  HI  and  obtain  average  correlations  for  various 
subgroups.  The  effect  of  source  iccation  on  the  autocorrelations  will  be 
examined  with  the  intention  of  gaining  information  about  propagation  mech¬ 
anisms.  Ensemble  average  and,  if  it  is  deemed  worthwhile,  subgroup 
average  autocorrelations  will  be  used  to  design  deconvolution  filters  for 
a  detailed  analysis  of  Ensembles  H  and  HI.  The  major  data  processing 
steps  are  summarized  in  the  block  diagram  in  Figure  H-17. 
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PROBABILISTIC  PROCESSING 


A.  INTRODUCTION 

Probabilistic  processing  is  a  method  of  processing  the 
output  data  of  an  array  of  seismometers  with  the  aim  of  detecting  earth¬ 
quake  or  underground  nuclear  blast  signals  in  the  presence  of  ambient 
seismic  noise.  This  method  is  based  on  the  assumption  that  the  array 
output  is  Gaussian  with  mean  zero  and  known  covariance  matrix  Q  or  0  , 
depending  on  the  absence  or  presence  of  signal. 

The  decision  regarding  presence  or  absence  of  the  signal 
is  therefore  made  equivalent  to  testing  the  hypothesis  that  the  observed 
data  are  from  a  Gaussian  population  with  covariance  matrix  Q  with  the 

simple  alternative  hypothesis  that  the  covariance  matrix  is  Q 

u  • 


The  univariate  situation  is:  a  single  value  x  is  observed 
from  a  Gaussian  population  with  mean  zero;  it  must  be  decided  whether  x 

2  2 

is  more  representative  of  a  Gaussian  population  with  variance  or  C  . 
The  two  density  functions  are 


V*) 

1 

exp  (-x2/20^) 

✓  3T  o. 

fp(x) 

1 

exp  (-x2/2o2) 

ii 

a| 
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ro 

III 


where  a  >  g  .  I£  the  original  likelihood  of  the  two  alternatives  were  equal, 

w  1 

it  is  more  likely  that  x  is  from  f  (x)  if  |xj  <  c  and  more  likely  that 

x  is  from  f  (x)  if  |x|  >  c.  So  the  method  of  processing  in  this  case  reduces 

to  squaring  the  observation  x  and  comparing  the  result  with  a  fixed  constant. 

If  a  vector  x  of  dimension  k  is  observed,  the  two  density 
functions  become 

Vx>  =  <2TT)”k/2  inj  ml/S  exp  (Jp*1  X/2) 

(2) 

f2(x)  -  (2rr)"k/2  |02l"1/2  exp  (-xT  O"1  x/2). 

These  equations  can  be  compared  to  find  the  set  of  vectors  x  for  which 
Mx)  >  f„(x).  Thus  if 

i  C. 

X^'O”1-  P”1)  x  >  c  =  In  (I  02 1 /!  I  )  (3) 


then  fg(x)  >  ^(x). 

T  - 1  - 1 

Therefore,  the  test  statistic  x  (0  ^  -  Q  )x  must  be 

computed  from  the  data.  The  detailed  calculation  of  the  test  statistic 
in  terms  of  array  data  is  given  in  the  following  section. 

A  more  rigorous  development  of  the  test  statistic  based 
on  Bayes  theorem  has  been  presented  in  a  preceding  special  report.  ^ 

B.  CALCULATION  OF  THE  TEST  STATISTIC 
1.  Notation 

The  data  from  an  array  cf  seismometers  can  be  represented 
by  the  matrix  X  where  the  rows  correspond  to  the  seismometers  and  the 
columns  to  che  time  sampled  data 


III -Z 


Thus,  x..  is  the  observation  from  seismometer  i  at  time  j. 
ij 

The  elements  of  matrix  X  can  be  "strung  out"  into  a  single 
vector  y  in  the  following  manner: 


,  T  T  Tv 

-  (x.,1  Xp,  .  .  .  Xj.; 


where  x.  =  (x,.,  x,, . ,  .  .  .  ,  x,  , . This  vector  vis  now  the  vector  of 

J  2j  Aj 

observations  from  which  the  test  statistic  (Equation  3)  must  be  computed. 
Since  the  mean  of  y  is  assumed  to  be  0,  the  covariance  matrix  of  y  is 

just  _ 

T 

0  =  yy 

y 

and  the  dimension  of  Q  is  MN  by  MN  since  there  are  MN  elements  in  the 

y 

vector  y.  Assuming  time  stationarity,  we  can  let 
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so  that 


0  = 
y 


where  each  submatrix  r.  is  an  M  by  M  matrix. 
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2.  Inverse  of  R 


so  that 


it  follows  that 


Now,  partition 


111-5 


N  .th  N 

Denote  by  s..  the  ij  M  by  M  submatrix  of  the  inverse  of  R  .  From 

Equations  (6)  and  (7),  the  following  equations  can  be  written: 


N 
s . . 
<J 


sn-!  +  rN'  (p.r1  rN'  T 

tj  N-i+1  VV  N-j 


'>  J 


j+1 

1,  .  .  N-i 


SN  .  =  s":1*  r1?  (P  r1  rN  T 

i+1,  j+i  ij  i+a  N  j+1 


i y  j  =  1;  •  •  • )  N-l 


(3) 

(9) 


Now,  by  subtracting  Equations  (8)  from  (9),  the  basic  iterative  formula 
of  Equation  (10)  is  obtained 


.N 

’i+1,  j+1 


=  s 


ij 


fN  (  rl  rN  T 
i+lv  N'  j+1 


-  rN'.  (p.  J-1  rN*  T  do) 
N-i+lV  N;  N-j+l  1  ' 


and  the  entire  inverse  matrix  can  be  generated  starting  with  the  first 

x  N 

M  columns  of  the  inverse  which  is  the  first  column  of  the  s..  's  available 

from  Equation  (5).  ^ 


3.  The  Quadratic  Processor 


Suppose  the  matrix  X  =  (x^,  .  .  .  ,  x^)  is  observed. 


Then  the  test  statistic 

P(v. 

i  ’  i  r  i 

can  be  evaluated  for  i  =  1,  ...»  T-N+l  where 


P(y.  )  =  y.CyT 


yj  -  (xT, 


C  -  P-'5'  -  P”"1 
L  1  2 


A  program  to  perform  chis  calculation  is  being  written.  In 
addition,  an  integrated,  squared  multichannel  Wiener  filter  output  will  be 
available  as  an  option  in  the  quadratic  processor  program  to  give  a  compar¬ 
ison  between  the  two  detection  methods. 
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4.  Multichannel  Vine  Series  Data  Generation 

It  is  desired  to  generate  multichannel  data  xT  =  (x^.,  x^,  •  ••>  ) 

N  pN 

with  the  covariance  matrix  R  .  Let  1  j  i  =  2,  .  .  .  ,  N  and  P  satisfy 

T  ” 

che  system  of  Equations  (5).  Suppose  y.  is  a  supply  of  vectors  satisfying 

y.=0 


- T  -  T 

y.  y.  -  i 


i  i 


Then,  x.  will  be  generated  by  the  following  system  of  equations: 


X,  =  €. 


X„  = 


1  *1  +  ‘2 

r2  *2  +  ^  *1  +  ‘3 


XN  =  r2  Vl  + 


*  +  ^  X1  +  €n 


(11) 


Xi  =  1 2  *i-i  +  ' 


‘  +  rN  Xi  -N  +  1  +  5i. 


The  method  of  determining  the  e.  is 


where  HH 


T 


€  =  Hy. 


The  reason  for  this  is  that 


€ .  must  satisfv 

i 


€  .  t  .  = 

i  I 


P 


N 


Determination  of  an  H  satisfying  the  above  is  aecornplished  by  requiring 
that  H  be  lower  triangular. 
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5.  Eigenvalues  (S'  igle  Channel) 

In  the  case  of  a  single  channel,  the  test  statistic  is 

P(x)  =  xV1  -  0:J)x  (12) 

-L  C. 

where  Q  ^  and  P  are  Toeplitz  matrices.  Let  the  observation  vector  x  be 
transformed  to  z  by  the  N  by  N  nonsingular  transformation  x  =  Sz  such  that 

P(x)  =  P ( z )  =  zVfP"1  -  P~2)  Sz  =  zT°z  (13) 

where  D  is  diagonal.  Then,  the  computation  of  the  quadratic  form  is 
reduced  to  filtering  operations  followed  by  a  summation  of  the  squared 

outputs.  An  especially  important  transformation,  diagonalizing  both  P  P 

simultaneously,  exists  since  A  =  P  and  ®  -  P0~^  are  positive  definite.  ^ 
This  transformation  is  found  by  solving  the  following  generalized  eigen¬ 
value  problem. 

The  solution  for  >  's  satisfying 

I A  -  XB1  =  0  (2U) 

is  called  the  generalized  eigenvalue  problem.  The  corresponding  set  of 
vectors  x  such  that 

Ax  =  XBx  (l^  ) 

are  the  generalized  eigenvectors.  The  matrix  M  whose  columns  consist  of 

T 

these  eigenvectors  normalized  so  that  x  Bx  =  1  then  can  be  shown  to 

satisfy  t 

M  AM  =  A 

MTBM  =  I  (.16  ) 

where  A  is  a  diagonal  matrix  with  elements  >  . 

Thus,  if  the  S  of  Equation  (13/  is  taken  to  be  M,  the 

diagonal  matrix  1)  satisfies 

D  -  I  -  A 

An  iterative  procedure  for  >  and  M  now  will  be  given.  Let 


f  ( X  )  -  | A  -  XB I 
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The  first  approximation  to  X  ,  the  smallest  eigenvalue,  is 


>  u  =  -f(o)/f  (0) 


(18) 


and  the  succeeding  iterative  values  are  found  by  the  chord  method.  Thus, 


f^Xl,  n-l|  Xl.  n-g  ~  f^Xit  n-2'  Xl,  n-1  , 


until 


is  satisfied. 


'In  ■  f(X.  )  -  w,  7^) 

a.,  n-i  l,  n-d 


X..  -  X  ^  <  X^  C  ;  r  - 

1,  r  1,  n-1  1>  n 

The  eigenvector  corresponding  to  X  is  computed. 


(19) 


20‘6  (20) 


It  is  desired  to  calculate 

f(X)/(X  - 


(21) 


in  order  to  continue  the  iterative  process.  Since  f(X)  is  assumed  to 
have  a  single  root  at  X  =  >  it  suffices  to  compute 


f  (>  )  - 


A  -  / 


l 


Z  8..  C.. 

i  j  i  ) 
tj 


(?.?.) 


where  C  is  the  matrix  of  cofactors  of  the  elements  of  A  ->  B.  The 
ij  1 

relation  (Equation  22)  is  shown  at  the  end  of  this  section. 


Now, 


(A.;  -  K  B..)  (C .  .  )  =  P 


1  Ij 


•J 


(23) 


since  >.  is  a  root  of  f('),  and  it  follows  that  each  column  (row)  of 
(C  )  is  proportional  to  the  unique  eigenvector  corresponding  to  X  . 

ij  1 

Thus,  ihe  rank  of  (C  )  is  one,  and  it  follows  that  the  matrix  (C.)  can 

ij  iJ 

be  written  in  the  form 


!cij)  =  b 


xx 


(?k) 


w 


here  x  is  the  eigenvector  corresponding  to  >  Thus, 


f  *  (X)L  ,  =  -I  B  C  =  -8xTBx  =  -8 

A=A]  IJ  IJ 


(25) 
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If  the  eigenvector  x  has  already  been  computed  , 
from 


r 

11 


S  can  be  determined 

(26) 


unless  C  j  ^  =  0. 
x  =  0,  so  that 


At  least  one  element,  say  j,  of  C^.  4  0,  for,  otherwise, 


Bx.x. 
1  J 


(27) 


would  determine  S  in  general. 

The  chord  method  again  is  used  to  find  the  smallest  0  of 
f(\)/(X->.),  the  first  approximation  of  X  being  X  and  the  second 

*  w  i  j  K  “  C* 

approximation  being  X^.  The  entire  set  of  solutions  X  and  eigenvectors 
x  are  solved  in  this  same  iterative  manner. 

The  relation  (Equation  22)  will  now  be  shown.  Write  out 

f(  > '  in  the  form 


Ali“XBli  A12“XB12  *  *  '  A1N"XB1N 


A2rXB2i 


,A2N_XB2N 


AN2"X3N . ANN"XBNN 


th 

If  the  >  in  the  ij  term  is  replaced  by  1 . .  where  X..  =  >  ,  then 

ij  ij 


fix',  =  z 

•  j 


Sf(XlI'  .  .  .  ,  XNN) 


ax. 


j 


Now, 


6f(Xn>  •  •  •  >  XW|J; 


N 


- - TT  (a.,  -  X  b  )  r  = 

ax  ax. .  ,  ,x  ik  ikik  ik  i  j  i j 

I  j  •  J  k=l 


and  Equation  (22)  follows. 
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PROGRESS  AND  FUTURE  GOALS 


n 

» 


Digital  programs  have  been  written  anti  checked  for 
computing  the  T  and  F*  matrices  and  the  quadra. ic  test  statistic.  The 

N 

program  for  the  inverse  matrix  R  has  been  written  but  not  checked  out. 
The  data  generation  program  has  not  been  written. 

It  is  expected  that  all  the  above  programs  will  be  completed 
by  50  May  1965  and  then  applied  to  1  eal  and  generated  data.  The  results 
will  be  presented  in  a  preliminary  report  by  15  July  1965. 

If  these  results  indicate  that  the  method  of  probabilistic 
processing  will  be  useful  in  the  problem  of  signal  detection,  a  program  to 
evaluate  the  eigenvalues  and  eigenvectors  of  the  relevant  covariance 
matrices  will  be  written  ai.d  the  data  used  to  reduce  the  evaluation  of 
the  test  statistic  to  filtering  and  squaring  over  a  gate.  This  more 
efficient  method  then  would  be  applied  to  a  set  of  real  data  and  the  results 
presented  in  the  December  1965  report. 
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SECTION  IV 


STUDY  OF  PARTIAL  ARRAYS  AT  CPO 
A.  INTRODUCTION 

Four  multichannel  filter  systems  for  CPO  were  designed 
using  the  full  array  and  three  partial  arrays.  These  systems  were  syn¬ 
thesized  using  measured  noise  correlations  and  an  infinite  velocity  signal 
model.  The  MCF  systems  were  evaluated  on  the  computer  to  compare 
the  ability  of  the  partial  arrays  to  suppress  noise  while  preserving  high- 
velocity  signals. 

Results  show  that  the  partial  arrays  are  almost  as  effective 
as  the  full  array  near  1  cps.  Below  approximately  0.7  cps,  seismometer 
gain  inequality  appears  to  be  effecting  the  output  of  the  systems.  At 
higher  frequencies,  the  K-space  aliasing  of  some  of  the  partial  arrays 
prevents  effective  noise  suppression. 

To  help  evaluate  the  effects  of  seismometer  gain  inequal¬ 
ities,  MCF  systems  will  be  designed  using  a  signal  model  incorporating 
statistical  gain  fluctuation.  This  will  be  the  subject  of  a  separate  report. 

b.  Multichannel  filtfr  systems  design 

The  filter  systems  were  synthesized  in  the  time  domain 
by  using  measured  noise  correlations  and  an  infinite-velocity  signal 
model.  In  all  cases,  desired  output  was  the  noise-free  signal  at  the 
center  seismometer.  All  filters  were  63  points  long. 

Noise  data  used  to  synthesize  the  systems  consisted  of  the 
averaged  sum  of  five  noise  samples  (A,  B,  E,  F,  I)  recorded  at  CPO  in 
1963.  These  noise  samples  had  been  prewhitened,  processed  through 
antialiasing  filters  and  resampled  to  a  sampling  interval  of  0.  07Z  sec. 

This  resulted  in  a  Nyquist  foldback  frequency  of  slightly  under  7  cps.  In 
this  report,  all  spectral  analysis  has  been  performed  on  the  range  of  0-7  cps. 
A  detailed  discussion  of  the  preparation  and  evaluation  of  these  noise  samples 

8 

has  been  published. 

Signal  correlations  used  in  the  time -domain  program  for 
designing  the  filters  were  generated  from  the  noise  power  spectrum  of 
of  the  center  seismometer.  The  power  spectrum  of  the  noise  used  in  the 
filter  design  program  was  scaled  up  by  4  to  give  the  signal  power  spectrum. 
Since  an  infinite -velocity  signal  model  was  used,  all  signal  auto-  and 
crosscorrclations  were  the  same. 
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No  statistical  gain  fluctuation  was  included  in  the  signal 
autocorrelations.  Results  indicate,  however,  that  seismometer  gain 
inequalities  are  affecting  the  MCF  systems.  This  will  be  discussed  in 
subsection  D. 

1.  CPO  IP  23 

This  system  consists  of  five  channels:  the  center  seis¬ 
mometer;  the  sum  of  three  seismometers  on  R2:  the  sum  of  three  seis¬ 
mometers  on  R3:  the  sum  of  six  seismometers  on  R4;  and  the  sum  of 
six  seismometers  on  R5  (a  total  of  19  seismometers).  See  Figure  IV-1. 

The  channels  which  are  sums  of  seismometer  outputs  are  not  averaged. 
Thus,  the  filter  designed  for  a  particular  ring  is  the  same  as  the  indi¬ 
vidual  filter  to  be  applied  to  each  seismometer  in  that  ring. 

2.  CPO  IP  24 

In  this  system,  the  center  seismometer,  the  sum  of  ring  3, 
the  sum  of  ring  4,  and  the  sum  of  ring  5  were  used  to  simulate  a  4-channel 
system  (16  seismometers). 

3.  CPO  IP  25 

For  this  system,  the  center  seismometer,  the  sum  of 
ring  4  and  the  sum  of  ring  5  were  used  to  simulate  a  3-channel  system 
(13  seismometers). 

4.  CPO  IP  26 

For  this  s/  stem,  the  center  seismometer  and  the  sum  of 
ring  5  were  used  to  simulate  a  2-channel  system  (7  seismometers). 

C.  EVALUATION  OF  THE  MULTICHANNEL  FILTER  SYSTEMS 

To  evaluate  and  compare  the  systems,  two  noise  samples 
and  an  artificial  infinite-velocity  signal  were  computer -processed.  The 
systems  were  analyzed  for  noise  suppression,  signal  preservation, 
signal/noise  improvement  and  wavenumbe r  response. 

The  noise  samples  used  in  the  evaluation  vere  CPO  noise 
samples  A  and  I.  They  were  recorded  over  a  month  apart  and  at  different 
times  of  the  day,  but  their  frequency  spectra,  sh  wn  in  Figure  IV-2,  are 
reasonably  similar.  These  spectra  and  all  othe^  power  soectra  were 
obtained  using  a  Tl-deveioped  spectral  progran  .  In  all  cases,  the  spectral 
plots  were  obtained  using  autocorrelations  with  time  lags  out  to  49. 

^  The  scales  of  all  power  spectra  in  this  report  are  related 

to  (rnu  '  /cps  at  1  cps.  In  this  report,  these  scales  are  correct  at  1  cps. 
The  p  ^whitening  filter  r  sponse  at  1  cps  has  been  taken  into  account. 

This  lias  not  been  the  case  with  previously  published  power  spectra. 


IV-2 


Figure  IV- 2.  Power  Density  Spectra  at  Center  Seismometer  —  Noise 
Samples  A  and  J. 
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Pre  nous  spectra  were  intended  to  be  related  to  ground  motion  at  any  fre¬ 
quency  with  the  aid  of  a  curve  incorporating  seismometer  PTA  response  and 
prewhitening  filter  response  (Figure  4-17,  Array  Research  Semiannual 
Technical  Report  No.  1).  Thus,  power  spectra  appealing  in  th;s  report 
will  be  shifted  up  3  db  from  those  previously  published. 

g 

From  the  above  cited  TI  report,  CPO  noise  is  shown  to 
have  the  following  characteristics: 

•  The  great  bulK  of  the  noise  power  is  concentrated 
below  4  cps 

•  The  noise  is  generally  quite  predictable  below 
2  cps,  indicating  tha'  CPO  noise  below  2  cps  is 
quite  coherent 

•  Noise  power  is  fairly  well  concentrated  in  K-space. 
Wavenumber  power  spectra  of  CPO  noise  at  ire- 
quencies  of  0.  5,  1,  2,  and  3  cps  are  shown  in 
Figures  IV-17,  IV-18,  IV-19,  IV -20  which  will  be 
discussed  under  subsection  F 


D.  MCF  SYSTEM'S  RESPONSES  TO  INFINITE- VELOCITY  SIGNAL 

The  amplitude  and  phase  response  of  IP  23,  24,  iP  25, 

and  IP  26  to  an  infinite-velocity  signal  are  shown  in  Figure  1V-3.  All 
systems  preserve  the  amplitude  of  an  infinite -velocity  signal  rather  well 
and  exhibit  near-zero  phase  re  'onse  wrhich  is  necessary  for  signal 
preservation. 

The  phase  and  amplitude  response  of  each  channel  of  the 
four  MCF  systems  are  shown  in  Figures  IV-4  through  IV-9.  The  phase 
response  of  some  of  these  individual  channels  shows  significant  deviations 
from  0  or  180°.  This  should  be  considered  in  detail  as  it  relates  to  the 
hypotheses  of  space-stationarity  and  seismometer  gain  equalization. 

Consider  a  MCF  system  designed  using  an  array  whose 
geometry  is  symmetrical  abou.  its  center.  The  wavenumber  response 
of  such  a  multichannel  filter  sy  stem  can  be  written  as 


where 
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Figure  IV -3 . 


Signa1  Responses  of  the  Multichannel  Filter  Systems 
Infinite  Velocity  Signal 


Fixing  f  and  taking  the  sum  over  pairs  of  symmetrically  placed  seis- 

— ♦  — *  — 4  — ♦ 

mometers  (call  these  locations  r  and  r  where  r,  =  r.,  =  -r  _),  the 

)l  j2  J  J 1  j2 

system  response  can  be  written  as 


+  H.0  e 
J2 


where 


-i2rr  |  k  |  |r  j  cos  0 

L,  h'  e  J 

j=i  ji 


i2n  k  r .  cos  9 


E  means  summation  by  symmetrically  placed  pairs 

j  =  l 

is  the  filter  applied  to  th*3  seismometer  placed 
nl  the  center  of  symmetry- -if  any 

•**4  •■4 

9  is  the  angle  between  k  and  r.^  (it  is  a  function  of 

r . ,  and  k  ) 

J1 

Thus,  if  the  filters  applied  to  symmetrically  spaced  seis¬ 


mometers  are  the  same,  the  system  response  becomes  2H.  cos 

|2:r  | k  |  |r.|  ^os  9j  +  H^.  This  can  be  considered  in  the  more  useful  form 

J 

E  a.  (k)  (r(H.)  +  I(H.) 

J=°  J  I  J  J  j 


where 


given  by 


a  is  real  and  not  a  function  of  frequency 

j 

th 

R|H.j  means  the  real  part  cl  the  j  filter 

th 

I|H.j  means  the  imaginary  part  of  the  j  filter 

The  MSE  of  ‘'he  MCF  system  at  any  fixed  frequency  is 


MSE  = 


=  ff  I  |l  E  a.  (k)  R  (H.J  +  [E  a.  (k)  I  (H  \1  j  S  (k)  dk  dk  + 

JJ  U  j  J  1  ,}J  h  3  1  J,J  I  x  y 


rr  r  -  2  i2) 

E  a.  (k)  R  +  E  a-  (k)  I  (H.\  (  N  (k)  dk  dk 

j  J  '  j'J  l  j  J  J,J  I  x  y 
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Figure  IV -4.  Amplitude  and  Phase  Responses  of  the  Filters  for  IP  23 
Channels  1  and  2 
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Figure  IV-8.  Amplitude  and  PI 
Channels  1 ,  2  an 
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Figure  IV-9.  Amplitude  and  Phase  Responses  of  the  Filters  for  IP  26 
Channels  1  and  2 


The  first  integral  gives  the  signal  distortion  part  of  the 
MSE,  with  S (k)  being  the  signal  power  distribution,  and  the  second  integral 

gives  the  output  noise  power,  N (k)  being  the  noise  power  distribution.  The 

2 


Ea  (k)  I  |H.j 


expression  containing  the  imaginary  parts  of  the  filters 

is  non-negative  as  are  S  (k)  and  N  (k) .  Thus,  the  MSE  is  minimized 
with  respect  to  I  when 


Since  the  a.(k 


's  are  linearly  independent  functions  of  k  ,  i.  e. 


a^(k)  =  2  cos  |2rrjkj  j  r .  |  cos  9 j 

the  above  can  be  true  over  a  k  region  only  if  I^H.j  =  0  for  all  j.  Therefore, 


given  an  MCF  system  whose  channels  were  obtained  from  sums  of  symmetri¬ 
cally  placed  seismometers,  it  would  be  expected  that  the  response  of  each 
filter  would  be  real  at  every  frequency. 


For  any  such  system  whose  filters  are  not  purely  real, 
the  hypothesis  that  the  signal  and  noise  are  plane  waves  is  not  satisfied. 
This  might  be  the  result  of  any  combination  of  the  following  factors: 

(1)  Signal  i"  not  space-stationary. 

(2)  Noise  is  not  space-stationary. 

(3)  Statistical  discrepancies  caused  by  finite  data 
analysis  make  the  sign'l  or  noise  appear  lonspace- 
stationary. 

(4)  Seismometer  gains  (or,  more  generally,  the 
frequency  responses)  of  symmetrically  placed 
pairs  are  unequal. 

CPO  systems  IP  25  and  IP  26  have  individual  channels 
which  are  the  sum  of  symmetrically  placed  seismometers  with  the  same 
filter  applied  to  each  seismometer.  (In  effect,  the  seismometer  outputs 
on  any  ring  are  actually  summed  before  filtering.)  If  the  noise  and  signal 
were  plane  waves,  the  above  would  insist  that  the  filter  response  of  each 
channel  be  real.  This  is  not  the  observed  response,  however.  (See 
Figures  IV -8,  IV-9.)  The  CPO  systems  in  this  report  were  designed 
using  artificial  plane  wave  signals.  Also,  previous  experience  indicates 
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that  the  effects  of  changes  due  to  statistical  variation  cf  the  data  should  be 
small.  Thus,  the  probable  causes  are  reduced  to  combinations  of  (2)  and 

(4). 

R5  of  IP  26  (Figure  IV-9)  shows  noticeable  deviation  from 
0  or  180°  phase  below  1/2  cps  and  in  the  range  3-7  cps.  At  the  higher 
frequencies,  local  crustal  changes  probably  tend  to  somewhat  destroy 
space  stationa rity.  Seismometer  gain  inequalities  probably  contribute 
to  the  nonreal  response  also.  Below  1  cps,  however,  it  is  more  reasonable 
to  infer  that  seismometer  gain  inequalities  are  the  dominant  factor.  Indeed, 
ambient  noise  generally  has  been  found  to  be  quite  space-stationary  at  low 
12 

frequencies . 

IP  25,  R4  and  R5  (Figure  IV-8),  exhibit  even  more  deviation 
from  0  or  180°  phase.  The  inference  is  the  same  as  for  IP  26.  This  system 
has  an  additional  ring  of  six  seismometers  as  compared  to  IP  26.  Seis¬ 
mometer  gain  inequalities  among  thesi  other  seismometers  could  tend  to 
make  the  individual  channels  of  this  MCF  system  more  phase-deviant.  The 
channels  might  have  to  become  more  phase-deviant  to  compensate  for  each 
other.  In  this  system  also,  nonreal  frequency  responses  of  the  channels 
are  marked  below  1  cps. 

CPO  MCF  systems  IP  23  and  IP  24  were  designed  on  arrays 
which  are  not  entirely  symmetric.  Thus,  a  directional  noise  field  could 
result  in  nonreal  frequency  response.  Frequency  responses  of  individual 
channels  of  these  systems  (Figures  IV -4  through  IV -7)  exhibit  such  general 
deviation  of  0°  or  180°  phase  response  that  seismometer  gain  inequality 
again  must  be  implied.  This  is  particularly  true  at  low  frequencies  where 
differences  in  seismometer  responses  on  a  given  ring  due  to  directional 
effects  would  be  minimized. 

Study  of  the  phase  responses  of  the  individual  channels  from 
the  various  systems  strongly  indicates  seismometer  gain  inequalities,  partic¬ 
ularly  below  1  cps.  K-space  responses  also  indicate  the  same  thing.  This 
will  be  discussed  under  subsection  F.  To  evaluate  the  effect  of  seismometer 
gain  inequalities  on  the  partial  arrays  further,  MCF  systems  incorporating 
statistical  gain  fluctuation  in  the  signal  model  will  be  synthesized  and  eval¬ 
uated.  This  will  be  the  subject  of  a  special  report  which  will  be  forth¬ 
coming  shortly. 

The  ability  of  the  four  systems  to  preserve  signals  of  less 
than  infinite  velocity  will  be  discussed  under  subsection  F. 
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E.  NOISE  SUPPRESSION  AND  SIGNAL-TO-NOISE  IMPROVEMENT 


Figure  IV -10  shows  time  functions  of  the  two  input  noise 
samples  and  the  filtered  output  of  each  MCF  system  for  each  noise  sample. 
These  MCF  systems'  output  traces  then  were  used  for  the  spectral  and 
spectral- ratio  analysis  presented  in  this  report.  The  5-ring  summed 
traces  of  each  noise  sample  show  the  noise  as  it  was  input  into  the  MCF 
systems.  The  traces  are  explained  in  detail  in  the  figure. 

IP  23,  24  and  25  appear  to  reject  the  noise  almost  equally 
well,  while  IP  26,  which  was  designed  using  only  the  center  seismometer 
and  the  outside  ring  of  six  does  a  noticeably  poorer  job. 

Figures  IV - 1 1  and  IV- 12  show  the  power  spectra  of  NSA 
processed  by  IP  23,  24,  25,  and  26.  Figure  IV- 13  shows  the  signal-to- 
noise  improvement  (signal  power  out/noise  power  out  div:‘  led  by  signal 
in/noise  power  in)  obtained  for  NSA  from  these  same  MCF  systems. 

Since  the  signal  responses  of  all  the  systems  are  so  nearly  uniform 
(Figure  IV-3),  S/N  improvement  ratios  are  essentially  a  measure  of 
noise  suppression. 

For  NSA,  the  S/N  improvements  obtained  from  IP  23,  24, 
and  25  show  only  moderate  differences  and  are  nearly  identical  up  to 
2-1/2  cps.  IP  26  shows  significantly  less  S/N  improvement  between 
1.  5  and  5  cps  but.  below  1.  5  cps,  gives  results  comparable  with  the 
other  systems. 

Figures  IV-14  and  IV-15  show  the  power  spectra  of  NSI 
processed  by  IP  23,  24,  25,  and  26.  Figure  IV-16  shows  the  S/N 
improvement  obtained  for  NSI  from  these  same  systems. 

Between  about  3/4  and  1-1/2  cps,  all  systems  give  similar 
S/N  improvement.  Below  3/4  cps,  S/N  improvement,  increases  as  the 
number  of  seismometers  increases.  Below  3/4  cps,  improvement 
obtained  from  IP  26  for  NSI  is  down  3-4  db  as  compared  with  IP  23.  This 
was  not  observed  for  NSA  (Figure  IV-13).  Indeed,  below  3/4  cps,  S/N 
improvement  is  generally  better  for  NSI  than  for  NSA  in  these  low  fre¬ 
quencies.  S/N  improvements  for  NSA  from  all  systems  are  comparable 
to  those  obtained  by  IP  26  for  NSI.  Any  combination  of  the  following  could 
account  for  this  discrepancy: 

•  Seismometer  frequency  responses  may  have  changed 
between  the  timt  NSA  was  recorded  and  NSI  was  recorded. 
They  were  taken  over  a  month  apart. 

•  The  seismic  sources  and  thus  the  K-space  distribution 

of  the  low-frequency  noise  are  different  iOr  NSA  and  NSI. 
These  two  noise  samples  show  a  somewhat  different  power 
spectra  below  3/4  cps. 
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Figure  IV-16.  Signal-To-Noise  Ratio  Improvements  Obtained 
Processing  "SI  and  an  Infinite  Velocity  Signal 


Above  1-1/2  cps,  the  MCF  systems  appear  to  operate 
on  NSI  in  a  manner  generally  similar  to  NSA.  IP  23,  24  and  25  give 
similar  S/N  improvement  to  2-1/2  cps.  above  2-1/2  cps,  the  differences 
among  these  three  systems  are  generally  2-4  db  but  range  to  over  10  db. 
Above  1-1/2  cps,  S/N  improvement  obtained  from  IP  26  is  generally 
down  5-10  db  as  compared  with  IP  23. 

From  a  comparison  of  S/N  improvements  obtained  from 
the  four  MCF  systems,  the  following  conclusions  can  be  reached* 

•  IP  25  (i3  seismometers)  generally  has  noise  rejection 
capability  nearly  equal  to  that  of  IP  23  (19  seismometers) 
below  2-1/2  cps. 

•  Below  1-1/2  cps,  IP  26  (seven  seismometers)  shows  noise 
rejection  which  is  generally  only  moderately  less  than  IP  23 
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•  Seismometer  gain  inequalities  may  be  giving  a  false  picture 
for  NSI  below  3/4  cps.  Greater  S/N  improvement  shown  by 
the  system  employing  more  seismometers  may  be  false  S/N 
improvement  obtained  by  using  the  seismometer  gain  in¬ 
equalities  to  suppress  noise.  The  artificial  signal  used  in 
obtaining  the  S/N  improvements  would  not  be  affected  by 
seismometer  gain  inequalities.  Thus,  such  a  system  would 
give  tetter  S/N  improvement  for  the  synthesized  situation 
than  for  an  actual  signal.  Below  1  cps,  seismometer  gain 
inequa  lities  are  indicated  by  the  wavenumber  responses  of 
the  systems  (subsection  F). 


Thus,  for  noise  suppression  in  the  frequency  range 
where  distant  mantle  P-wave  power  is  concentrated, 

•  Both  IP  25  and  26  are  generally  satisfactory. 

•  IP  23  is  not  a  significantly  better  system  than 
IP  25. 

•  IP  23  is  only  moderately  superior  to  IP  26. 

F.  K-SPACE  RESPONSE  OF  SYSTEMS 

Wavenumber  responses  of  the  four  systems  were  computed 
at  0.  5,  1.  0,  2.  0,  and  3.  0  cps.  A  mathematical  development  of  this  tech¬ 
nique  can  be  found  in  John  Burg's  "Three-Dimensional  Filtering  with  an 
Array  of  Seismometers ,  "  Geophysics,  October  1964,  p.  699-700. 

Results  are  presented  as  contour  maps  in  wavenumber  space  (Figures 
IV-17,  IV-18,  IV-19,  and  IV-20).  Contours  below  approximately  -24  db 
were  removed  to  make  a  more  presentable  map. 

K-space  responses  of  the  four  systems  at  a  particular 
frequency  have  been  combined  into  one  figure  with  tne  wavenumber  power 
spectra  of  the  noise  at  that  frequency.  K-space  power  spectra  of  the  noise 
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are  reproductions  of  those  previously  published. 

Figure  IV-17  shows  the  wavenumber  responses  of  IP  23, 

24,  25,  and  26  and  the  noise  power  spectra  at  0.  5  cps.  Comparison  of 
wavenumber  responses  and  noise  power  spectra  indicates  seismometer 
gain  inequalities  at  this  frequency.  The  noise  power  is  concentrated  in  an 
annulus  representing  Velocities  around  2.5  k.n/sec.  IP  23  attenuates  most 
strongly  plane  waves  with  a  velocity  of  about  2.  5  km/sec.  IP  24  attenuates 
most  strongly  plane  waves  with  velocity  approximately  2.6  km/sec. 
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IP  25  attenuates  most  str^rgly  plane  waves  with  a  velocity  of  about  3  km /sec, 
and  IP  26  attenuates  most  strongly  plane  waves  with  a  velocity  of  about 
3,5  km/sec.  The  systems  show  markedly  different  K-space  responses. 

This  is  characteristic  of  MCF  systems  which  have  seismometer  gain 
inequalities.  The  K-space  plots  show  the  MCF  systems  responses  as 
they  would  appear  if  all  seismometers  were  gained  equally.  The  actual  MCF 
systems,  which  incorporate  seis  gain  inequalities,  probably  have  similar 
wavenumber  responses  at  0.5  cps.  A  detailed  study  of  seismometer  gain 
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inequalizal  c.i  effects  has  been  published. 

The  noise  power  at  0.  5  cps  falls  off  rapidly  on  the  annulus 
representing  horizontal  velocities  of  less  than  2  km/ sec.  This  strongly 
affects  the  responses  of  the  systems.  Indeed,  IP  26,  which  suppresses 
noise  almost  as  well  as  the  other  systems  at  this  frequency,  actually 
gains  up  the  power  of  plane  wave  with  these  velocities. 


At  0.  5  cps,  all  systems  preserve  power  in  plane  waves 
with  velocities  greater  than  12  km/sec.  In  no  case  is  such  a  wave 


attenuated  by  as  much  as  3  db. 

Figure  IV-18  shows  wavenumber  responses  of  IP  23,  24, 

25,  and  26  and  the  noise  power  spectra  in  the  K-plane  at  1  cps.  Compar- 
ison  of  noise  spectra  with  the  various  system  responses  gives  good  general 
agreement  with  the  S/IJ  improvement  of  the  multichannel  filter  systems 

at  1  cps. 


On  the  disk  representing  velocity  greater  than  3  km/sec, 
ail  the  MCF  systems  have  similar  responses.  This  corresponds  to 
velocities  where  the  noise  power  is  generally  concentrated.  As  wou 
be  expected,  all  systems  give  about  the  same  S/N  improvement  (10  db) 
at  1  cps  with  a  difference  of  only  about  1  db. 


Aliasing  and  pseudo-aliasing  show  up  in  the  MCF  response 
maps,  but  their  effect  is  outside  the  range  of  appreciable  noise  power 
concentration.  The  only  true  alias  shows  on  the  map  of  IP  2b. 


Ail  MCF  systems  preserve  plane  wave  energy  with  hori 
zontal  velocities  greater  than  12  km/sec.  On  the  circle  representing 
these  velocities,  power  is  down  less  than  3  db  for  all  systems. 


Figure  IV-  19  shows  wavenumber  responses  of  IP  23,  24, 

25  and  26  and  the  noise  power  spectrum  in  K-space  at  2  cps.  Compari 
s™  of  system  responses  to  noise  iistribution  gives  generally  good  agree- 
ment  with  the  observed  S/N  improvement  of  the  MCP  systems  at  ^  cp  . 


annulus 

annulus, 

K-spare 


The  noise  power  at  2  cps  is  strongly  concentrated  in  an 
representing  velocities  between  3.  5  and  5  km/sec.  In  this  same 
IP  23,  24  and  25  exhibit  attenuation  of  lb  db  or  more.  In  genera 
response  of  each  of  these  systems  is  verv  similar.  At  2  cps, 


u  -  3  3 


. . . . . . . . . . 


channels  2  and  3  of  IP  23  and  channel  3  of  IP  24  are  contributing  relatively 
little  to  the  output  of  ch”  MCF  systems.  At  2  cps,  all  these  systems  show 
equal  S/N  improvement  within  1  db. 

IP  26  shows  space -aliasing  which  results  in  passing  much 
of  the  noise  energy  at  2  cps.  As  this  would  indicate,  S/N  improvement 
at  2  cps  is  about  12  db  less  for  IP  26  than  for  the  other  systems.. 

IP 23,  24  and  25  preserve  plane  wave  power  for  velocities 
above  12  km/sec.  These  systems  show  attenuation  of  4  db  or  less  on 
this  range  of  velocities.  IP  26,  however,  attenuates  plane  waves  with 
a  velocity  of  12  km /sec  approximately  10  db. 

Figure  IV -20  shows  wavenumber  responses  of  IP  23,  24, 

25,  and  26  and  the  spectrum  of  the  noise  power  at  3  cps.  Comparison  of 
system  responses  with  noise  power  spectra  again  gives  generally  good 
agreement  with  the  observed  S/N  improvement  ratios  at  3  cps. 

The  noise  power  is  generally  concentrated  in  an  annulus 
representing  velocities  around  4  km/sec.  Total  noise  power  at  3  cps 
is  down  5-10  db  from  the  frequencies  previously  considered. 

In  the  stronger  noise  band,  the  responses  of  IP  23,  24 
and  25  are  quite  variable.  IP  23  attenuates  the  noise  somewhat  better 
than  IP  24  which,  in  turn,  is  a  better  system  than  IP  25.  True  aliasing 
occurs  in  IP  25  just  on  the  ^dae  of  the  band  of  noise  power  concentration. 

At  3  cps  the  S/N  improvement  differences  among  these  three  systems 
are  about  5  db,  with  IP  23  being  best  and  IP  25  poorest  on  both  processed 
noise  samples. 

Close-spaced  aliasing  dominates  the  IP  26  wavenumber 
response  at  3  cps.  This  system  would  have  relatively  poor  ability  to 
attenuate  noise  in  any  reasonable  velocity  band.  As  this  would  indicate, 
the  S/N  improvement  of  IP  26  is  down  about  10  db  from  that  of  IP  23 
at  3  cps. 

At  3  cps,  none  of  the  MCF  systems  satisfactorily  preserves 
the  energy  of  12-km/sec  plane  waves.  All  systems  show  attenuation  of  more 
than  6  db  in  this  velocity  range.  This  bad  tendency  becomes  more  marked 
as  the  number  of  seismometers  is  reduced.  IP  26  attenuates  high-velocity 
plane  waves  (V>  12  km/sec)  as  much  as  18  db. 

Analysis  of  wavenumber  responses  of  the  various  systems 
taken  with  S/N  improvement  indicates  that 
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•  Seismometer  gain  inequalities  appear  to  be 
affecting  MCF  systems  significantly  below  1  cps 

•  Below  about  1-1/2  cps,  IP  26  (seven  seismometers) 

is  almost  as  good  a  system  as  IP  23  (19  seismometers); 
i.  e.  ,  below  1.5  cps,  IP  26  preserves  signals  with 
apparent  horizontal  velocities  >  12  km/ sec  and  suppresses 
noise  almost  as  well  as  IP  23. 

•  Below  2-1/2  cps,  IP  25  (13  seismometers)  is  very 
nearly  the  equal  of  IP  23  (19  seismometers).  Again, 
the  desirable  system  preserves  signal  with  propa¬ 
gation  velocities  >  12  km/ sec  and  suppresses  noise. 

Thus,  for  evaluation  of  distant  P-wave  signals,  IP  25  is 
essentially  as  good  as  IP  23,  and  IP  26  is  nearly  as  good  as  IP  23.  It  is 
apparent  that  R2  and  R3  contribute  little  to  the  ability  of  the  CPO  array 
to  enhance  distant  mantle  P-wave  signals.  R4  is  of  questionable  value. 

Seis  gain  inequalities  may  be  affecting  the  comparisons  below  3/4  cps, 
so  the  value  of  R4  has  not  been  definitely  determined.  MCF  systems 
which  will  be  designed  incorporating  statistiral  gain  fluctuations  should 
help  evaluate  the  usefulness  of  R4. 
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SECTION  V 


INITIAL  PROCESSING  AND  ANALYSIS 
OF  THE  LONG  TFO  NOISE  SAMPLE 


This  section  contains  information  concerning  the  recording 
and  preliminar  y  nr oces sing  and  analysis  of  a  20-min  ambient  seismic  noise 
sample  from  the  ring  array,  the  large  cross  array  and  a  3-component  seis¬ 
mometer  at  Tonto  Forest  Seismological  Observatory.  Also  presented  is  an 
analysis  of  calibration  data  for  the  large  cross  array  recorded  on  the  same 
day  as  the  long  noise  sample.  Section  VI  contains  a  more  detailed  analysis 
of  this  noise  sample. 


A.  RECORDING  AND  PRELIMINARY  PROCESSING  AND  ANALYSIS 

A  long  sample  of  ambient  seismic  noise,  slightly  less  than 
20  min  long,  was  edited  from  a  DFS  recording  taken  at  Tonto  Forest  Seis¬ 
mological  Observatory  on  December  20,  1963.  During  this  time  period, 
average  wind  velocity  was  1  to  5  mph,  and  the  temperature  was  20°F.  This 
noise  sample  was  free  of  recording  errors  and  appeared  to  be  free  of  non- 
seismic  noise  such  as  spiking.  The  sampling  period  of  the  field  recording 
was  24  msec,  which  gives  50,  000  points  in  20  min.  Because  of  limitations 
on  input  computer  capacity,  this  long  noise  sample  was  divided  into  10 
samples  labeled  C,  D,  E,  F,  G,  H,  I,  J,  K,  and  L,  respectively,  which 
were  2  min  each  in  duration  with  a  2  sec  overlap  in  time  for  every  two 
adjacent  noise  samples.  Thus,  each  of  the  10  noise  samples  contained 
5000  points  with  an  average  overlap  of  83-1/3  points  in  the  adjacent  noise 
samples.  The  GCT  start  time  for  noise  sample  C  was  05:55:00,  andthe 
GCT  stop  time  for  noise  sample  L  was  06:14:42. 

Figure  V-l  shows  the  geometry  of  Tonto  Forest  Seismological 
Observatory.  Table  V-l  shows  the  trace  position  of  seismometers  and  rings 
on  the  field  reel.  The  first  five  traces  contain  the  ring-averaged  summations. 
Ring  1  is  seismometer  Z16,  and  the  remaining  four  rings  arc  the  averaged 
sums  of  the  seismometers  in  their  respective  rings.  Thus,  ring  2  consists 
of  seismometer  Z10,  Zll,  Z17,  Z22,  Z?  1 ,  and  Z15;  ring  3  of  seismometers 
Z5,  Z12,  Z23,  Z27,  Z20,  and  Z9;  ring  4  of  seismometers  Z4,  Z6,  Z18, 

Z28,  Z26,  and  Z14;  and  ring  5  of  seismometers  Zl,  Z2,  Z7,  Z13,  Z24, 

Z29,  Z31,  Z30,  Z25,  Z19,  Z8,  and  Z3.  Traces  6  through  16  are  the  11 
seismometers  in  the  NS  arm  of  the  cross  array,  with  the  seismometer s 
going  in  order  from  south  to  north.  Traces  17  through  26  are  the  seismom¬ 
eters  in  the  EW  arm,  proceeding  east  to  west  except  that  the  center  seis¬ 
mometer  is  not  repeated  and  that  seismometers  70  and  25  are  interchanged. 
Trace  27  contains  vertical  seismometer  Zl.  Traces  28  and  29  contain  the 
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figure  V-l.  Arrays  at  Tonto  Forest  Scismological  Observatory  (TFO) 
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Table  V-l 


TRACE  IDENTIFICATION,  MOTOR  CONSTANTS 
AND  GROUND  DISPLACEMENTS 


Trace 

Seismometer 

G 

Ground  Displacement 

No, 

Identification 

(newtons /amp) 

(mu  0-P) 

1 

Z16 

2 

R2 

3 

R3 

4 

R4 

5 

R5 

6 

Z63 

0.425 

17.92 

7 

Z62 

0.  410 

17.29 

8 

Z61 

0.430 

18.13 

o 

Z60 

0.420 

17.71 

10 

Z9 

0.415 

17.50 

li 

Z21 

0.405 

17.08 

12 

Z31 

0.435 

18.35 

13 

Z71 

0.415 

17.50 

14 

Z72 

0.445 

18.77 

15 

Z73 

0.425 

17.92 

16 

Z74 

0.420 

17.21 

17 

Z67 

0.415 

17.50 

18 

Z66 

0.415 

17.  50 

19 

Z65 

0.395 

16.  66 

20 

Z64 

0.  350 

14.  76 

21 

Z13 

0.415 

17.50 

22 

Z17 

0.445 

18.77 

23 

Z70 

0.410 

17.29 

24 

Z68 

0.420 

17.71 

25 

Z69 

0.405 

17.08 

26 

Z25 

0.425 

17.92 

27 

Z1 

28 

Z36  E-W 

29 

Z37  N-S 

30 

Cal.  T.  B. 

31 

GC.T 
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cast-west  and  north-south  horizontal  seismometers,  respectively.  These 
two  seismometers  are  located  in  the  vault  with  seismometer  Zl.  Trace  30 
which  is  used  for  recording  the  calibration  signal  during  station  calibration 
was  blank  during  this  recording,  and  trace  31  contained  the  GCT  timing 
pulses . 


The  tirst  27  traces,  i.e.  the  five  rings,  the  linear  cross 
array  and  seismometer  Zi,  of  each  of  the  10  noise  samples  were  filtered 
wit*  a  symmetric  43-point  prewhitening  and  antialiasing  filter  and  then 
decimated  by  3  to  yield  filtered  72-msec  data.  Traces  28  and  29  which 
contain  the  output  of  the  ease-west  and  north-south  horizontal  seismometers, 
respectively,  were  not  filtered  but  were  decimated.  The  square  wave  func¬ 
tion,  which  was  placed  in  trace  30  during  the  running  of  the  Edit  II  program 
and  may  be  used  to  generate  tinting  lines  for  playback  purposes,  also  was 
decimated  by  3.  The  Edit  II  information  on  these  10  noise  samples  has  been 

2 

presented  in  a  previous  report.  The  frequency  response  oj  the  prewhiten  - 
filter  is  in  Figure  V-2.  In  applying  this  43-point  filter  to  the  10  noise  sam¬ 
ples,  only  the  compl  ?tely  filtered  outputs  were  desired.  Thus,  the  outputs 
were  31  points  shorter  on  each  end.  This  was  the  reason  for  having  overlap 
in  the  10  noise  samples.  To  join  the  10  noise  samples  back  into  one  long 


Figure  V-2.  Frequency  Response  of  Prewhitening  and 
Antialiasing  Filter 


4 


sample  after  decimation  by  3,  special  care  had  to  be  taken  in  where  to  start 
the  decimation  of  each  noise  sample  to  insure  that  the  long  noise  sample  did 
contain  every  third  point  in  time  after  joining  the  decimated  samples.  When 
the  10  filtered  and  decimated  noise  samples  were  appended,  there  were 
lb,407  points  with  a  frequency  range  from  0  to  6-17/18  cps. 

Power  density  spectra  were  comptited  from  the  long  noise 
sample  for  each  of  the  21  seismometers  in  the  linear  cross  array.  Figures 
V-3  through  V-6  iihow  these  spectra  which  were  obtained  from  autocorrela¬ 
tions  with  lags  cut  to  t  =  ±61.  Also,  the  spectral  ratios  relative  to  the  cen¬ 
ter  seismometer,  seismometer  Z21,  were  computed  for  each  of  the  other 
20  seismometers.  These  results  are  pictured  in  Figures  V-7  through  V-10. 

Two  channel  coherence  estimates  also  were  computed  at  this 
time  for  seismometers  Z63  and  Z64,  Z65  and  Z66,  and  Z68  and  Z69.  These 
coherences,  shown  in  Figures  V-ll  and  V-12,  were  found  from  correlation 
functions  with  lags  out  to  t  =  ±  124.  In  this  set  of  illustrations,  the  notation 
63  y  ^4  refers  to  the  coherence  between  seismometers  Z63  and  Z64. 

The  prewhitened  long  noise  sample  was  bandpass -filtered  and 
again  resampled  by  3  to  yield  216-msec  data  with  a  frequency  range  of  0  to 
2-17/51  cps.  At  this  point,  the  data  also  were  scaled  by  a  factor  of  4.  Fig¬ 
ure  V-13  shows  the  frequency  response  of  this  long  decimation-by-3  filter 
over  the  range  of  0  to  6-17/18  cps  and  also  of  0  to  2-17/54  cps. 

A  6-point  deconvolution  filter  designed  from  seismometer 
Z21  with  10  percent  white  noise  was  applied  to  the  216-msec  data.  The  fre¬ 
quency  response  of  this  filte  .  j  shown  in  Figure  V-14,  and  the  combined 
response  of  the  three  filters  which  have  been  applied  to  the  long  noise  sam¬ 
ple  are  shown  in  Figure  V-15.  For  completeness.  Figure  V-16  shows  the 
response  of  the  J-M  PTA  system  used  at  TFO. 

Power  density  spectra  with  maximum  lag  of  T  =  ±61  for  seis¬ 
mometers  Z62,  Z21,  Z73,  Z65,  and  Z70  are  shown  in  Figure  V-17.  These 
spectra  were  taken  after  the  216-msec  data  were  deconvolved. 

B.  CALIBRATION  ANALYSIS  OF  THE  LARGE  CROSS  ARRAY 

A  calibration  rur.  at  1  cps  on  the  linear  cross  array  was 
recorded  oi  the-  same  date  as  the  preceding  noise  recording.  The  GCT 
start  time  for  the  calibration  was  04:15:32,  and  the  GCT  stop  time  was 
04:16:00  for  all  21  seismometers.  Taking  into  account  the  motor  constants 
of  the  seismometers,  these  calibration  data  v/ere  analyzed  to  give  absolute 
scales  for  the  power  density  spectra  of  the  earth  motion  at  1  cps. 
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Figure  V-6.  Power  Spectra  of  Seismometers  Z13,  Zl7,  Z70, 
Z68,  Z69,  and  Z25 
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igure  V-7.  Power  Spectra  Ratios  of  Z63,  Z62  .  Z6l,  Z60,  and  Z9 
Relative  to  Seismometer  Z21 
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Figure  V-8.  Power  Spectra  Ratios  of  Z31,  Z7i,  Z72,  Z73,  and  Z74 
Relative  to  Seismometer  Z21 
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Figure  V-9.  Power  Spectra  Ratios  of  Z67,  Z66,  Z65,  Z64,  and  Z13 
Relative  to  Seismometer  Z21 
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Figure  V-10.  Power  Spectra  Ratios  of  Z17,  Z70,  Z68,  Z69»  and  Z25 
Relative  to  Seismometer  Z21 
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Figure  V-l  1 . 


Coherences  Between  Seismometers  Z66  and  Z65 
and  Between  Seismometers  Z68  and  Z69 


On  occasion,  the  motor  constant  of  calibration  coil  G  in  the 

J-jv1  seismometers  at  TFO  is  checked.  The  motor  constant  is  determined 

in  the  following  manner.  Several  weight  lifts  are  taken  with  a  weight  whose 

mass  is  W  grams,  and  the  trace  deflection  X1Ir  is  recorded  in  millimeters 
l  w 

for  each  trial.  Then  d-c  pulses  I  in  amperes  are  applied  to  the  calibration 

coil  and  the  current  adjusted  until  the  trace  deflection  due  to  current  X.  is 

i 

within  10  percent  of  X^  in  amplitude.  The  corresponding  X.  and  I  then  are 
used  in  the  following  formula  to  calculate  G: 


FREQUENCY  (cps) 

Figure  V-12.  Coherence  Between  Seismometers  Z63  and  Z74 


newtons  /ampere 


Table  V-l  also  shows  the  motor  constants  of  the  seismometers  in  the  linear 
cross  array.  These  motor  constants  were  taken  from  a  data  sheet  dated 
December  16,  i  96 3 . 
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Figure  V-14.  Frequency  Resoonse  of  the  Deconvolution  Filter 


The  formula  used  for  computing  earth  motion  in  millimicrons 
zero-to-peak  is 

G  x  I  x  106 

Earth  Motion  in  mp.  0-P  -  T ~~z 

4  r  f  M 

where  G  is  the  motor  constant  in  newtons /ampere,  I  is  the  input  current  in 
milliamps  0-P,  f  is  the  frequency  in  cps,  and  M  i«*  the  mass  of  the  seismom¬ 
eter  in  kilograms.  The  mass  of  the  Johnson-Matheson  seismometer  was 
18.02  kg,  frequency  was  1  cps  and  input  current  was  0.03  ma  0-P.  The 
computed  equivalent  earth  motions  also  are  listed  in  Table  V-l. 

Ground  displacement  for  seismometer  Z63  was  17,  92  mp 
zero-to-peak.  At  the  time  of  the  calibration,  the  DFp  contained  a  -12  db 
attenuator  card  but,  during  the  noise  recording,  no  attenuation  was  used. 
Therefore,  effective  displacement  can  be  considered  to  be  4.48  mp  zero-to- 
peak. 
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Figure  V-16.  Response  of  the  J-M  Data  PTA  Seismometer  System 
Used  at  TFO 
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The  transfer  function  of  the  seismometer  to  magnetic  tape 

B 

system  at  i  cps  is  given  by  K  -  —  where  A  is  the  zero-to-peak  equivalent 

A 

earth  motion  in  millimicrons  at  1  cps  and  B  is  the  zero-to-peak  amplitude 
of  the  recorded  numerical  cosine  function  on  the  tape. 

The  number  B  is  found  from  a  computer  program  which 


•  Determines  frequency  f  of  the  calibration  signal 
(recorded  on  trace  30)  by  counting  the  zero-crossing 
of  the  sine  wave 

•  Calculates  the  Fourier  transform  of  the  seismom¬ 
eter  traces  g  (t)  at  that  frequency,  i.e.  , 


N  -  1 

C  (f)  =  At  V]  g  ( j A t )  (cos  2tt f j A t  -  i  sin  2nfjAt) 
n  '  n 

j  =  0 


B  then  is  given  by 


2lC(f) 

NAt 


where  A  t  is  the  sampling  period. 


This  program  also  calculates  the  d-c  level  of  the  traces  and  the  amplitude 

and  phase  of  C  (f)  relative  to  trace  11  which  is  the  center  seismometer  in 
n 


the  cross  array.  The  output  of  this  program  is  shown  in  Table  V-2  and 
explained  in  Table  V-3. 


(null  2 

The  0  db  relative  to  1  - *—  at  1  cps  is  equal  to  1  J  log,  AK 

cps  1 0 


=  20  log^gK.  Hence,  for  seismometer  Z63,B  =  ryy 


2(36554] 


(1 1 38 )(0.  024) 


and,  thus. 


201°S10 

All  the  0-db  values  for  the  21  seismometers  were  computed  according  to 
the  preceding  process. 

Table  V-4  gives  the  0-db  ground  motion  value  for  the  center 
seismometer  in  the  cross  array  and  the  relative  0-db  values  for  the  other 
seismometers.  Also  listed  is  10  times  the  logarithm  of  the  avciagc  square 
value  of  the  tape  numbers  for  the  center  seismometer  during  various  pro¬ 
cessing  stages  and  the  relative  db  values  for  the  other  traces.  Variations 
in  the  relative  db  value  for  a  trace  after  the  various  processing  stages  are 
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Table  V-2 


CALIBRATION  ANALYSIS  OUTPUT  FOR  TFO  CROSS  ARRAY 


(a)  Calibration  on  Linear  Cross  Array  at  TFO 

(b)  001747 
001016 

(c)  000001 

(d)  001138 

(e)  000024 

(f)  000027 

(g)  000030 

(h)  .  98828E+0 

(i)  .  20667E+3 


(j) 

(k) 

(1) 

(m) 

(n) 

000006 

000032 

.  36554E+5 

. 10153E+1 

-.  22031E+1 

000007 

000032 

.  34223E+5 

. 95071E+0 

.  05245E+0 

000008 

000028 

.  36263E+5 

.  1  0073E+1 

-.  24094E+1 

000009 

000024 

.  36392E+5 

.  10109E+1 

-.  37767E+0 

000010 

000024 

.  34517E+5 

. 95886E+0 

-.  34411E+1 

000011 

000016 

.  35997E+5 

. 1 0000E+1 

-.  00000E+0 

000012 

000008 

.  36834E+5 

.  1 023 lE-t  l 

-.  36230E+1 

000013 

000003 

.  35809E+5 

.  99477E+0 

-.  13254E+1 

000014 

-000004 

. 36161E+5 

.  10046F+1 

-.  30773E+1 

000015 

-000004 

. 37787E+5 

.  1 0498E+1 

-.  26682E+1 

000016 

-000009 

.  37855E+5 

.  10515E+1 

27556E+1 

000017 

000007 

.  38848E+5 

. 10791E+1 

65394E40 

000018 

000091 

.  38261E+5 

. 10629E+1 

-.  20666E+1 

000019 

000012 

.  36289E+5 

. 1 0080E+1 

-.  21401E+1 

000020 

000014 

.  35863E+5 

. 9q626E49 

,  261  57E+1 

000021 

00001  0 

.  35962E+5 

.  "9902E+0 

.  182  54E+1 

000022 

000015 

.  37569E+5 

.  1043  7E+1 

. 34235E+1 

000023 

000027 

.  37588E+5 

.  10441E+1 

.  10162E+2 

000024 

000021 

.  36351E+5 

.  10097E+1 

.  1 9582E+1 

000025 

000031 

.  38456E+5 

.  10683E  +  1 

.  35180E+1 

000026 

000020 

.  36766E+5 

.  10212E+1 

. J4516E+1 
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Table  V-3 


EXPLANATION  OF  TABLE  V-2 


(a)  Title,  assigned  on  input 

(b)  Input  event  number 
Output  event  number 

(c)  Block  number  of  first  zero  crossing  of  calibration  trace 

(d)  Number  of  blocks  chosen  as  Fourier  transform  gate 

(e)  Sample  rate 

(f  )  Number  of  cycles  in  gate 

(g)  Trace  on  which  calibration  appears 

(h)  Computed  frequency  of  the  calibration  signal 

(i  )  Phase  angle  of  reference  trace  in  degrees 
(j  )  Trace  number 

(k)  D-C  level 

(l )  Absolute  value  of  Fourier  transform 

(m)  Relative  amplitude 

(n)  Relative  phase  angle  in  degrees 


due  to  differences  in  its  power  spectrum  relative  to  the  spectrum  of  the  cen¬ 
ter  seismometer  (trace  11).  These  variations  are  small.  However,  in  com¬ 
paring  the  relative  total  powers  and  the  relative  0- db  levels,  there  are 
appreciable  differences  on  some  traces.  In  particular,  the  calibration  on 
seismometer  Z9  (trace  10)  says  that  its  noise  power  is  more  than  3  db  less 
than  Z21  (trace  11).  These  results  show  at  least  that  the  calibration  mea¬ 
surements  are  a  poor  way  of  normalizing  ambient  seismic  noise. 
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Table  V-4 


O-DB  GR  DUND  MOTION  VALUES  AND  TOTAL  POWER  OF  TRACES 
AFTER  VARIOUS  PROCESSING  STAGES  FOR  TFO  CROSS  ARRAY 


0-db  Trace  Power  in  db  Trace  Power  in  db  Trace  Power  in  db  after 

Ground  after  Prewhitening  before  Prewhitening  Repeat  Decimation  by  3, 

Motion  and  Decimation  by  3  and  Decimation  by  3  Gain  by  4  and  Deconvolution 


6 

-0.3 

-0.4 

0.2 

0.  1 

7 

-0.6 

-1.3 

-0.4 

-0.4 

8 

-0.5 

-0.8 

-1.4 

-1.2 

9 

-0.2 

-0.5 

-1.0 

-1.0 

10 

-0.6 

-3.8 

-4.5 

-4.3 

11 

55.  8 

61.3 

43.7 

65.5 

12 

-0.4 

-0.4 

-1.0 

-1.2 

13 

-0.3 

-0.6 

-0.2 

0.0 

14 

-0.8 

-0.7 

-0.2 

0.6 

15 

0.  0 

-0.3 

-0.6 

-0.4 

16 

+0. 1 

-0.  1 

-0.3 

-0.1 

17 

+0.  4 

0.2 

0.5 

0.6 

18 

+0.  3 

0.3 

0.3 

0.6 

19 

+0.3 

-1.7 

-1.8 

-1.3 

20 

+1.2 

-2.  0 

-1.9 

-1.4 

21 

-0.2 

-1.9 

-1.8 

-1.3 

22 

-0.5 

-1.  5 

-1.6 

-1.2 

23 

+0.  3 

-2.0 

-1.8 

-1.7 

24 

-0.2 

-1.6 

-2.0 

-2.0 

25 

+0.  6 

0.1 

-0.6 

-0.4 

26 

-0.2 

0.  2 

0.2 

0.4 
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SECTION  VI 


FURTHER  ANALYSIS  OF  TFO  AMBIENT  NOISE 

A.  SUMMARY 

This  section  presents  the  complete  set  of  all  auto-  and  cross¬ 
correlation  functions  from  the  21  seismometers  in  the  large  cross  array  at 
TFO  as  found  from  the  long  noise  sample  discussed  in  Section  V.  Also  pre¬ 
sented  is  the  complete  set  of  10  coherences  between  the  two  horizontal  seis¬ 
mometers  and  three  vertical  seismometers.  The  analysis  of  this  data  shows 
that  the  4  to  6  see  energy  is  coming  from  N60°  E  *vith  a  velocity  of  about 
3.4  km/sec.  The  ocher  energy  below  1,2  cps  is  mantle  P-wave  noise,  which 
is  definitely  non-isotropic.  Above  1.2  cps,  the  energy  is  random,  seismo¬ 
meter  to  seismometei  ,  except  for  some  highly  coherent  lines  in  the  spectra. 
The  horizontal  seismometer  in  line  with  the  4  to  6  sec  energy  shows  a  high 
coherence  with  the  vertical  seismometers  in  that  frequency  range,  but  the 
transverse  seismometer  has  very  small  coherence  with  all  other  seismo¬ 
meters.  This  is  further  proof  that  this  energy  is  Rayleigh  motion.  The 
presence  of  Love  wave  energy  on  the  transverse  seismometer  can  be  inferred 
from  its  lack  of  coherence  with  the  other  seismometers. 


B.  PROCESSING  OF  THE  AUTO-  AND  CROSSCORRELATION  FUNCTIONS 

All  possible  auto-  and  crosscorrelation  functions  between  the 
21  seismometers  in  the  large  cross  array  at  TFO  were  computed  from  the 
long  ambient  noise  sample  discussed  in  Section  V.  The  correlations  were 
computed  from  the  noise  traces  after  they  had  been  resampled  down  to  216  msec 
data  and  had  been  deconvolved.  Thus  the  Nyquist  fold -over  frequency,  W, 
was  2.  315  cps  and  the  power  spectra  were  almost  white  except  for  the  lack 
of  energy  below  0.  1  cps  and  the  presence  of  line  spectra.  The  correlations 
were  all  computed  out  to  ±6i  lags  or  13.  176  sec.  In  addition,  the  auto  -  and 
crosscorrelations  were  appropriately  scaled  to  make  the  zero  lags  on  all 
autocorrelations  equal.  (Some  noise  traces  were  low  in  power;  sec  Section  V.) 

Four  symmetrical,  41 -point,  digital  bandpass  filters  were 

W  W  W 

computed  by  Fourier  transformation  of  the  frequency  bands,  0  to  ,  T"  to  , 
WWW  t  o  3 

—  to  —  ,  and  —  to  W  (0  to  0.  386  cps,  U.  386  to  0.  772  cps,  0.  772  to  1. 157  cps, 

and  1.  157  to  2.  315  cps).  The  frequency  responses  of  these  4  bandpass  filters 
are  shown  in  Figure  VI-1.  Each  of  the  231  correlation  functions  were  filtered 
by  each  bandpass  filter.  Only  fully  filtered  outputs  were  calculated  and  thus 
the  resulting  correlation  functions  were  20  points  shorter  on  each  end;  i.  e.  , 
t  =  -t 4 1  lags  or  ±8.856  sec.  The  filtered  correlation  functions  were  not 
normalized  among  themselves,  but  the  results  from  the  first  three  filters  were 
scaled  up  by  6  and  the  results  from  the  last  filter  were  scaled  up  by  2  so  that 
the  resulting  autocorrelations  would  be  comparable  in  magnitude  to  the  unfiltered 
autocorrelations. 
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Figure  VI-1.  Frequency  Responses  of  the  Bandpass  Filters  Applied 
to  the  Coi  relation  Functions 

Figures  VI-2  through  VI -6  show  all  auto  -  and  crosscorrelations 
from  just  the  north -south  line  in  the  cross  array.  The  autocorrelations  lie 
along  the  main  diagonal  of  the  figures  and  go  in  order  from  left  to  right,  south 
to  north  as  can  be  seen  from  the  inset  <T  che  TFO  geometry.  Figure  VI-2 
shows  the  unfiltered  correlations;  Figure  VI-3,  the  0  to  0.386  cps  correlations 
Figure  VI-4,  the  0.  386  to  0.772  cps  correlations;  Figure  VI-5,  the  0.772  to 
l.  157  cps  correlations;  Figure  VI-6,  the  1.157  cps  to  2.315  cps  correlations. 
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Figure  VI -2.  Correlation  Functions  for  the  NS  Line  Full  Frequency  Band  of  0  to  2.  315  Cps 


Correlation  Functions  (X6)  for  the  NS  Line  Frequency  Band  of  0  to  0.  386  Cps. 


Correlation  Functions  (X6)  for  the  NS  Line  Frequency 
to  0.  77  2  Cps 
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Figure  VI-6.  Correlation  Functions  (X2)  for  the  NS  Line  Frequency  Band  of 
1.  157  Cps  to  2.  315  Cps 


The  crosscorrelations  in  these  figures  and  in  the  other  correlation  displays 
lie  at  the  appropriate  grid  locations.  For  example,  in  Figure  VI -2  the  upper 
right  hand  correlation  is  Z63XZ74.  The  definition  of  Z63XZ74 

1  N 

Z63XZ74(t)  =-  E  g63(n'  t)  g74(nAt-T) 
i  =  1 

is  used  for  the  crosscorrelation  functions. 

Figures  VI-7  through  VI-11  are  a  similar  set  for  the  east-west 
line  of  the  cross  array,  where  left  to  right  along  the  diagonal  is  east  to  west 
along  the  line.  Note  that  the  center  seismometer,  Z21,  is  repeated  in  the 
figures.  Figures  VI-12  through  VI-16  are  a  similarly  filtered  set  for  the 
correlations  between  the  two  lines  of  seismometers.  Again  the  center  seis¬ 
mometer  correlations  are  repeated  in  these  figures  for  spatial  consistency 
in  the  figures. 

C.  ANALYSIS  OF  THE  CORRELATION  DISPLAYS 

An  immediately  noticeable  property  of  the  correlation  functions 
is  their  space  stationarity.  All  autocorrelation  functions  are  quite  similar, 
reflecting  the  similarity  of  their  power  spectra.  However,  more  interesting 
is  the  similarity  of  the  crosscorrelations  which  lie  on  lines  parallel  to  the 
main  diagonal  in  Figures  VI-2  through  VI-11.  The  seismometers  in  the  two 
lines  are  approximately  equal  spaced  (about  l  km  apart)  and  thus,  if  the 
ambient  noise  field  at  TFO  is  space  station  '"y,  then  the  observed  similarity 
should  exist.  In  the  case  of  Figures  VI-12  to  VI-16,  space  stationarity 
requires  that  the  crosscorrelations  symmetrically  placed  with  respect  to 
Z21XZ21  should  oe  time  reverses  of  each  other,  which  is  approximately 
observed.  Considering  the  quality  of  the  data,  the  following  analysis  of  the 
wavenumber  structure  of  ambient  noise  at  TFO  is  on  a  firm  foundation. 

In  looking  at  the  low  frequency,  0  to  0.  386  cps  uand,  it  is 
evident  that  there  is  a  strong  directional  4  to  6  sec  period  wave  going  across 
the  array.  Figure  VI -3  shows  a  moveout  of  about  1,6  sec  across  the  10  km 
length  of  the  NS  line  with  the  energy  hitting  seismometer  Z74  before  Z63. 
Figure  VI -8  shows  a  smaller  moveout  across  the  length  of  the  EW  line  with 
the  energy  hitting  the  east  end  before  the  west  end.  The  energy  must  be 
strongly  directional  because  of  the  high  coherence  across  the  array.  In 
addition  to  this  low  frequency  energy,  there  is  evidence  of  higher  frequency 
energy  which  has  little  or  no  moveout;  i.  e.  ,  mantle  P -wave  energy  Because 
of  the  interference  of  these  two  types  of  energy,  Figure  VI-13  is  most  useful 
in  determining  the  direction  of  arrival  of  the  4  to  6  sec  energy.  By  looking 
at  the  correlations,  say,  in  the  last  column  of  Figure  VI-13,  the  Z60XZ20 
crosscorreiation  function  is  the  most  symmetric,  indicating  that  the  low 
frequency  wave  front  arrives  at  these  two  seismometers  at  the  same  time. 
From  this  pair  and  other  pairs  with  good  symmetry,  the  wave  motion  is 
consistently  found  to  lie  along  the  line  N60°E.  Estimating  the  moveout  between 
a  pair  of  seismometers  lying  along  this  line;  i.e.,  Z74  and  Z68  or  Z63  and 
Z13  (both  have  a  separation  of  5.8  km  and  about  1.7  sec  moveout)  a  velocity  of 
about  3.4  km/sec  coming  from  N60E  can  be  estimated  for  the  4  to  6  sec 
period  energy.  This  result  is  in  agreement  with  previous  work;  see  Section 
IT  of  reference  1  .  The  velocity  is  very  reasonable  for  Rayleigh  wave  motion. 
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Correlation  Functions  for  the  EW  Line  Full  Frequency  Band  of 
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Figure  VT-9.  Correlation  Functions  (X6)  for  the  EW  Line  Frequency  Band  of 
0.  386  to  0.7 7  2  Cps 
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Figure  VI-10.  Correlation  Functions  (X6)  for  the  EW  Line  Frequency  Band  of  0.77 
to  1.157  Cps 
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Figure  VI-11.  Correlation  Functions  (X2)  for  the  EW  Fine  Frequency  Band  of 

1.  157  to  2.  315  Cps 


lation  Funct.ons  Between  the  Two  Lines  Full  Frequency  Band  of 
.315  Cps 
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Figure  VI-14.  Correlation  Functions  (X6)  Between  the  Two  Lines  Frequency  Band  of 
0.  386  to  0.  772  Cps 
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Figure  VI-16.  Correlation  Functions  (X2)  Between  the  Two  Lines  Frequency  Band  of 


The  frequency  bands  of  0.  386  to  0.77Z  cps  and  0.77Z  to  1 .  1  57  cps 
appear  to  be  mainly  mantle  P-wave  energy  since  the  correlations  have  very 
little  rnoveout.  A  striking  phenomena  is  the  vanishing  and  reappearance  with 
phase  reversal  of  the  crosscorrclations  with  increase  in  separation  distance. 
This  property  indicates  a  range  of  velocities  for  the  mantle  P-wave  energy. 

In  particular,  if  the  energy  were  isotropically  distributed  from  **■  to  a  lowest 
velocity  V  (like  the  solid  disk  models  used  for  signals  in  designing  multi¬ 
channel  filter  systems),  the  null  distance,  X,  would  occur  when 


or 


JL(ZnX^)  =  0 


J  (1.  Z197  tt)  =  0 


V  = 


ZXf 
1.  Z197 


(See  reference  1Z.) 


(1) 


From  Figure  VI-9,  the  correlation  between  Z65  and  Z70  is  small, 
giving  a  null  distance  of  about  8.  5  km.  Using  a  center  frequency  of  .  58  for 
the  0.  386  to  0.77Z  cps  band,  the  value  of  V  given  by  Equation  (1/  is  about 
8  km/sec.  Along  the  NS  line  in  this  band,  Figure  VI-4,  the  null  distance  is 
the  full  length  of  the  line  or  about  10.7  km.  This  gives  a  V  of  about  10  km/sec. 
In  the  frequency  band  of  0.77Z  to  1 .  1  57  cps,  Figure  VI-10  nulls  at  Z17XZ70 
or  a  distance  of  5.4  km.  Using  a  center  frequency  of  0.97  cps,  the  velocity 
V  is  given  as  8.6  km/sec.  Figure  VI-5  nulls  at  Z61XZ70,  giving  a  V  of 
13.6  km/sec.  Thus,  along  the  north-south  arm  the  correlations  seem  to 
agree  with  a  solid  disk  model  with  a  minimum  velocity  of  about  1Z  km/sec, 
but  along  the  east -west  line  the  minimum  velocity  is  about  8  km/sec.  These 
results  would  fit  a  model  of  a  solid  elipse  with  the  longer  axis  along  the 
east-west  line.  Furthermore,  there  is  a  slight  rnoveout  on  the  correlation 
functions  from  one  end  to  the  other  end  on  both  lines  of  the  cross  array.  The 
moveouts  indicate  more  energy  is  coming  from  the  NNW,  From  these 
observations,  it  is  clear  that  the  distribution  of  mantle  P-wave  noise  is  not 
isotropic. 

Ir.  the  frequency  band  above  1 .  I  57  cps,  there  is  little  coherent 
energy  and  no  evidence  of  space  stationar ity .  The  only  coherent  energy  seems 
to  be  due  to  lines  in  the  spectra  of  the  seismometers  and  this  coherence  is  not 
space  stationary.  This  indicates  that  there  is  no  particular  velocity  structure 
to  the  ambient  noise  at  TFO  above  1  cps  which  is  in  agreement  with  previous 

results . 
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D.  PREPARATION  OF  THE  HORIZONTAL- VERTICAL  COHERE*  'CE  DATA 


The  noise  samples  from  the  EW  and  NS  horizontal  seismometers, 
which  are  located  in  the  vault  with  the  vertical  seismometer  Zl,  were  not  fil  - 
tered  by  the  prewhitening,  antialiasing  filter  used  on  the  rest  of  the  traces 
before  they  were  resampled  down  to  72  msec  data.  Refer  to  Section  V.  A 
digital  filter  was  designed  so  that  it  had  the  same  frequency  response  as 
the  prewhitening  filter  in  the  range  of  0  to  2-17/54  cps  and  applied  to  the 
216  msec  horizontal  traces.  Thus,  the  resulting  horizontal  traces  were 
corrected  so  that  the  total  digital  filtering  applied  to  them  agreed  with  that 
applied  to  the  vertical  seismometers. 

Because  the  4  to  6  sec  energy  was  coming  from  N60°E,  the 
two  horizontal  seismometers  were  combined  to  produce  (1)  a  rotated  EW 
seismometer,  REW,  which  would  be  in  line  with  this  microseismic  noise; 
and  (2)  a  rotated  NS  seismometer,  RNS,  which  would  be  transverse  to  this 
energy.  Thus,  the  linear  combinations  used  were 

REW  =  0.  5  NS  +  0.  366  EW 


and 

RNS  =  0.  866  NS  -  0.  5  EW  (2) 

Three  vertical  seismometers,  Zl,  Z13,  and  Z25,  were 
selected  to  be  analyzed  in  connection  with  the  rotated  horizontals.  This 
set  of  verticals  gave  one  seismometer  located  at  the  two  horizontals,  one  in 
line  with  the  low  frequency  directional  energy  at  a  distance  of  2.  2  km  and 
one  transverse  to  the  energy  at  a  distance  of  3  km.  After  all  filtering  had 
been  completed,  a  small  portion  of  the  noise  samples  for  Z13,  Z25,  Zl, 

EW,  NS,  RNS  and  REW  is  shown  in  Figure  VI-17. 

Figure  VI-18  shows  the  auto -power  spectra  of  the  three 
vertical  seismometers  and  the  two  rotated  horizontal  seismometers.  Two 
auto -power  spectra  estimates  for  each  seismometer  were  obtained  from  the 
two  channel  coherence  program,  one  using  correlations  with  lags  of  ±61 
(13.  176  sec)  and  the  other  with  lags  out  to  ±124  (26.784  sec).  Calibration 
data  for  the  vault  seismometers  was  lacking,  so  that  there  is  no  absolute 
scales  given  for  these  spectra.  Figures  VI-19  through  VI-23  show  the 
complete  set  of  10  coherences  between  the  five  seismorr.  2rs  for  the  two 
different  correlation  lag  values. 

E.  ANALYSIS  OF  THE  HORIZONTAL- VERTICAL  COHERENCES  AND  POWER 
SPECTRA 

In  looking  at  Figures  VI -18  thru  VI-23,  the  high  coherence 
between  the  inline  horizontal  with  the  three  vertical  seismometers  for  the 
4  to  6  sec  energy  is  striking.  The  transverse  horizontal  has  no  coherence 
with  the  verticals  or  with  the  inline’  horizontal  except  for  line  spectra  in  the 
higher  frequency  regions. 
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Portion  of  the  Noise  Samples  for  Z13,  Z25,  Zi,  EW,  NS,  RNS  and  REW 


UNCALIBRATED  POWER  SPECTRA  (db) 
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Figure  VI-21.  Coherences  For  ZIXREW  and  Z1XRNS 
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This  is  a  clear  indication  that  the  4  to  6  sec  energy  is  Rayleigh  wave  energy 
and  is  coming  from  about  N60°E.  The  difference  in  the  spectra  of  the  two 
horizontal  seismometers  (Figure  VI -18)  in  the  frequency  range  of  0.  1  to 
0.5  is  open  to  discussion.  However,  the  lack  of  coherence  of  the  energy  on 
the  transverse  horizontal  v/ith  t.  other  seismometers  indicates  that  it  is 
Love  wave  energy  coming  from  the  same  direction  as  the  Rayleigh  wave 
energy. 

These  results  indicate  that  an  array  of  3-component  seis¬ 
mometers  could  be  used  to  separate  teleseismic  signals  from  Rayleigh  and 
Love  wave  ambient  noise. 
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